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Fourier  analysis  of  the  approximation  power  of 
principal  shift-invariant  spaces 
Carl  de  Boor  &  Amos  Ron 


1.  Introduction 


Spaces  spanned  by  finitely  or  countably  many  translates  of  one  or  several  basic  functions  play 
an  important  role  in  spline  theory,  radial  basis  function  theory,  sampling  theory  and  wavelet  theory. 
Spline  theory  stresses  the  case  when  the  basic  functions  are  compactly  supported,  while  sampling 
theory  singles  out  the  case  when  the  spectrum  (i.e.,  the  support  of  the  Fourier  transform)  of  the 
basic  functions  is  compact.  In  the  radial  basis  function  theory,  neither  of  these  is  assumed,  and 
instead,  the  computational  simplicity  as  well  as  the  positive  definiteness  (i.e.,  the  positivity  of 
the  Fourier  transform)  of  the  basic  functions  is  preferred.  Finally,  wavelet  theory  focuses  on  the 
interrelation  between  the  initial  space  and  its  dyadic  dilates.  In  all  these  areas,  the  underlying 
space  s  is  meant  for  approximation  or  decomposition  of  functions,  and  thus,  the  determination  of 
its  approximation  properties  is  of  basic  significance.  : . 

The  present  literature  is  mainly  concerned  with  a  space  s  which  is  the  algebraic  or  topological 
span  of  the  integer  translates  of  one  basic  function  ip.  More  precisely,  we  hold  a  collection  {iph}hei 
of  complex- valued  measurable  functions  defined  on  fit*1,  where  I  is  either  the  open  interval  (0 . .  h0), 
or  a  discrete  subset  of  such  an  interval  (e.g.,  {2"n  :  n  6  IN}).  For  each  h,  we  look  at  all  linear 
combinations  Ylaeh'Z*  ^('~a)a(a)>  f°r  which  this  sum  converges  in  a  certain  sense,  and  denote  by 
Sh  the  space  of  all  limit  functions  obtained  in  this  way.  Roughly  speaking,  we  call  Sh  the  span  of  the 
h7Zd -translates  of  ip  hi  and  this  is  an  exact  description  of  $h  in  case  iph  is  of  compact  support,  a  case 
in  which  the  sum  £Q6/l2Zd  ^/i(-  ~  <*)<*(<*)  *s  locally  finite,  and  hence  arbitrary  linear  combinations 
are  allowed  in  this  sum.  Approximation  properties  are  primarily  studied  via  approximation  orders: 
for  the  given  scale  {sh}hi  one  examines  the  rate  of  decay  (as  h  -*■  0)  of  dist(/,s/,),  where  /  varies 
over  some  space  of  admissible  functions,  which  must  contain  all  test  functions  in  2?(IRd)  (namely, 
all  C00(IR<i)  compactly  supported  functions),  and  the  distance  dist (f,Sh)  between  /  and  $h  is 
measured  in  some  norm,  usually  a  p-norm  (1  <  p  <  oo),  or  a  weighted  p-norm.  We  say  that  the 
approximation  order  of  the  scale  {s/,}/,  is  k  (or  0(hk)),  for  some  positive  (usually  integer)  k , 
if,  for  every  admissible  /,  dist(/,s/l)  =  0(hk ),  with  a  constant  that  depends  on  /  (and  clearly  not 
on  h),  while,  for  some  admissible  /,  dist (f,Sh)  £  o(hk). 

Although  a  discussion  of  the  above  model  can  already  be  found  in  Schoenberg’s  work  [S] 
(for  unhariate  functions),  the  first  comprehensive  analysis  of  approximation  orders  was  carried  out 
about  twenty  years  ago  primarily  by  people  from  the  finite  element  group,  the  best  known  reference 
for  which  is  [SF].  Strang  and  Fix  considered  the  “compactly  supported  scaling  case”,  namely,  when 
ipi  is  compactly  supported  and  iph  is  its  /i-dilate  (i.e.,  iph  =  ipi(-/h)),  and  showed  (for  the  2- 
norm)  that  approximation  orders  are  characterized  by  the  polynomials  in  sj.  Some  modifications 
and  improvements  of  these  results  (known  these  days  as  “The  Strang-Fix  Conditions”)  can  be 
found  in  [DM2],  [BJ]  and  [R2].  However  [DR],  the  polynomials  in  spline  spaces  are  unrelated  to 
approximation  orders  if  the  {iph}h  are  not  the  dilates  of  one  function.  Discussion  of  approximation 
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orders  for  compactly  supported  piecewise-exponentials  {4>h}h  can  be  found  in  [DR],  [BR]  and  [LJ]. 
We  know  of  no  study  of  approximation  orders  for  general  compactly  supported  {il>h}h- 

In  the  study  of  the  above  problem,  one  usually  considers  separately  the  questions  of  lower 
bounds  and  upper  bounds  on  the  approximation  order  (and  hopes  of  course  to  match  them).  The 
standard  approach  to  lower  bounds  is  via  the  quasi-interpolation  argument:  first,  a  space  H  C  Pl^s/, 
is  identified,  and  then  the  local  approximation  properties  of  H  are  converted  to  approximation 
orders  of  {$/,}/,  with  the  aid  of  local  linear  operators  (=quasi-interpolants)  whose  restriction  to 
H  is  the  identity.  The  space  H  consists  of  polynomials  in  the  scaling  case,  and  of  exponential- 
polynomials  in  the  piecewise-exponential  case,  but  need  not  to  be  so  in  general  (cf.  [BAR]).  Further, 
the  condition  H  C  s/,,  all  h,  is  convenient,  but  not  essential,  as  the  quasi-interpolation  argument 
of  [R2]  shows.  For  earlier  constructions  of  quasi-interpolants  see,  e.g.,  [SF]  and  [BF].  An  updated 
discussion,  together  with  a  partial  bibliography,  can  be  found  in  [B2]  and  [BR]. 

In  contrast  to  lower  bounds,  there  does  not  seem  to  exist  a  standard  approach  to  the  upper 
bound  question.  We  already  mentioned  [SF]  and  [BJ],  and  we  add  [LC],  [JL]  and  [HL],  where 
weaker  forms  of  approximation  orders  (“local”,  “controlled-local”)  are  characterized,  under  the 
assumption  that  the  basic  functions  are  either  compactly  supported  or  maintain  a  high  order  0  at 
oo  (where  “high”  is  defined  relative  to  the  desired  approximation  order,  and  several  basic  functions 
are  allowed  in  each  h-layer).  However,  all  these  results  are  confined  to  the  scaling  case,  and  further, 
the  fast  decay  at  oo  that  is  required  from  the  basic  functions  excludes  various  functions  of  interest. 
Sharp  upper  bounds  on  the  approximation  order  of  polynomial  box  spline  spaces  and  exponential 
box  spline  spaces  (integer  direction  case)  were  derived  in  [BH]  and  [LJ]  respectively,  based  on  the 
local  structure  of  the  spline  space,  which  in  general  is  a  rare  possibility  (see  the  box  spline  section 
in  this  paper).  Optimal  schemes  for  approximating  bounded  exponentials  in  the  non-scaling  (still, 
compactly  supported)  case  were  introduced  in  [R?].  These  results  will  be  presented  in  the  sequel, 
since  they  form  the  starting  point  for  the  upper  bound  analysis  here. 

We  introduce  and  analyze  in  this  paper  a  new  approach  for  the  determination  of  the  approxi- 
mation*orders  of  the  scale  {s/Ja.  In  this  approach,  only  modest  decay  rates  are  required  of  the  basic 
function  iph  (e.g.,  some  maximal  function  ifif  should  be  integrable),  and  the  questions  of  upper 
bounds  and  lower  bounds  are  attacked  almost  simultaneously,  so  that,  for  all  special  cases  studied 
here,  they  match  each  other  and  the  approximation  order  is  determined.  Using  Fourier  analysis 
methods,  we  further  need  not  restrict  our  attention  to  integral  approximation  orders.  On  the  other 
hand,  for  the  lower  bound  part,  we  place  some  smoothness  conditions  on  the  basic  functions,  which 
are  met  in  all  examples  we  know  from  the  radial  basis  function  theory,  but  exclude  splines  of  low 
smoothness,  so  that  we  have  here  the  usual  smoothness-localization  trade-off.  This  approach  makes 
no  use  of  quasi-interpolation  arguments;  in  particular,  polynomial  or  exponential  reproduction  is 
not  required.  In  addition,  the  approximation  scheme  is  constructive  enough  for  the  determination 
of  realistic  estimates  for  the  constant  which  is  hidden  in  the  0(hk)  expression. 

In  spite  of  the  generality  of  the  results  here,  we  are  able  to  apply  them  directly  to  obtain  upper 
and  matching  lower  bounds  for  the  case  when  the  basic  function  is  a(n  exponential)  box  spline  with 
rational  directions.  We  believe  that  none  of  the  methods  now  in  the  literature  could  provide  either 
bounds.  We  show  the  important  fact  that  many  of  the  lower  bounds  known  for  radial  basis  (and 
related)  functions  underestimate  the  correct  approximation  order,  and  explain  this  phenomenon. 
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Finally,  we  show  that  the  use  of  basic  mollifiers  for  the  basic  function  (e.g.,  the  Gaussian  kernel) 
leads,  if  properly  used,  to  infinite  approximation  orders. 

As  mentioned,  lower  bounds  on  the  approximation  order  were  derived  previously  with  the 
aid  of  quasi-interpolants,  and  the  difficulty  we  observed  in  the  implementation  of  this  method 
encouraged  us  to  start  the  work  reported  here.  While  the  quasi-interpolation  argument  is  an 
extremely  useful  and  powerful  tool  in  the  compactly  supported  scaling  case,  its  application  in 
other  known  situations  is  complicated.  For  example,  for  piecewise-exponentials,  the  space  H  of 
exponential-polynomials  in  C)h3h.  might  be  hard  to  determine,  its  local  approximation  properties 
might  be  even  harder  to  analyze  (cf.  [DR]),  and  the  lower  bounds  attained  in  this  way  might 
underestimate  the  true  approximation  order  (all  these  three  are  valid  difficulties  in  the  exponential 
box  spline/rational  direction  case).  But  the  major  drawback  of  the  quasi-interpolation  argument 
appears  in  the  area  of  radial  basis  functions  (cf.  [P]  and  the  references  therein).  In  almost  all 
examples  there,  iph  is  the  h-dilate  of  ip\,  hence  one  expects  to  use  polynomial  reproduction  in 
the  quasi-interpolation  argument.  Still,  if  the  function  ipi  does  not  decay  fast  enough,  standard 
polynomial  reproduction  arguments  (namely,  Poisson’s  summation  formula)  do  not  apply.  Further, 
even  if  all  desired  polynomials  are  shown  to  be  reproduced,  more  subtle  ‘^formation  on  the  rates  of 
decay  of  ip j  is  required,  [DJLR],  [Bul-3].  At  the  outset  of  our  present  study,  we  tried  to  apply  to 
these  cases  the  quasi-interpolation  argument  from  [R2],  which  involves  only  bounded  exponentials, 
but  found  that,  although  the  polynomial  reproduction  argument  can  be  circumvented  in  this  way, 
no  better  approximation  orders  are  obtained. 

The  crux  of  all  the  analysis  here  is  the  linkage  between  the  Fourier  transform  and  Fourier 
series  via  the  periodization  argument,  and  which  is  best  expressed  by  Poisson’s  summation  for¬ 
mula.  Starting  with  [S],  this  tool  has  always  been  the  chief  Fourier  analysis  argument  for  polyno¬ 
mial/exponential  reproduction.  The  results  of  this  work  show  that  the  periodization  argument  is 
not  only  an  important  technical  tool,  but  is  at  the  center  of  the  approximation  order  analysis:  the 
rearrangement  of  the  error  into  Fourier  series  allow-:  us  to  distinguish  between  terms  that  can  be 
reduced  by  an  optimal  selection  of  the  approximant,  and  terms  that  can  be  small  only  because  of 
the  good  approximation  properties  of  the  spaces  {s/,}/,. 

We  have  chosen  in  this  paper  to  focus  on  the  L <„  case,  namely,  measure  the  error  in  the 
oo-norm,  primarily  since  this  substantially  simplifies  the  analysis  of  upper  bounds  (by  making  the 
exponential  functions  admissible  for  approximation).  On  the  other  hand,  this  norm  is  probably  one 
of  the  harder  choices  in  the  lower  bound  analysis  (certainly  when  compared  to  the  2-r.orm):  Since 
the  approximation  is  performed  entirely  in  the  Fourier  transform  domain,  we  needed  to  bound  the 
Sobolev  (or  potential)  norm  of  the  function  to  be  approximated  in  terms  of  its  Fourier  transform, 
and  thus  our  notion  of  “admissibility”  falls  short  of  the  usual  Sobolev  space.  Further,  as  we  already 
mentioned  before,  our  lower  bound  conditions  exclude  basic  functions  of  low  smoothness,  and  this 


is  again  related  to  the  choice  of  the  norm:  the  error  in  the  approximation  scheme  can  be  written 
and  analyzed  in  terms  of  certain  Fourier  multipliers,  whose  Fourier  transform  is  explicitly  known. 
However,  to  obtain  sharp  results  with  the  aid  of  these  multipliers  requires,  because  of  the  use  of 
the  oo-norm,  information  about  the  behaviour  of  the  multiplier  in  the  original  domain,  which,  as  a 
rule,  is  not  easily  accessible. 

Throughout  the  paper,  C  stands  for  the  unit  cube  [-1/2 ..  1  /2]n,  and  Br,  for  the  /^(IR^-ball 
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of  radius  t)  centered  at  the  origin.  We  use  the  notation  e$,  6  G  IR4*,  for  the  complex  exponential 


and  denote  by  <j> *'  the  semi-discrete  convolution 

<(>*':/»  £  ^(*  “  <*)/(«)> 
o£2 l* 

where  /  is  any  function  defined  (at  least)  on  7Ld.  The  Fourier  transform  of  the  summable  function 
/  is  defined  by 

f(0)  :=  /  e-e(t)f(t)  dt, 

Jr* 

and  is  extended  by  duality  to  all  distributions  in  V(TR.d).  We  also  make  use  of  the  discrete 
Fourier  transform  (or  symbol)  /  of  the  function  /  (of  polynomial  growth),  defined  as 

/:=  £  «-./(»). 
aOL* 


Note  that 

(1.1)  fW  =  (/*'e.)(0)  =  (e../*'l)(0) 

in  case  /j  4  6  li(7Ld).  We  denote  by  II  the  ring  of  all  polynomials  in  d  variables,  and  IIn  is  the 
subspace  of  polynomials  of  degree  at  most  n.  Also,  II-n  ;=  IIn_i. 

As  a  rule,  a,  /?  are  generic  points  of  7Zd,  2ir7Ld,  respectively,  and  0 ,  w  are  generic  points  of  the 
Fourier  transform  domain.  Also,  the  default  norm  is  ||  •  II  •—  II  •  ||oo>  while,  for  x  G  IR<<,  |x|p  is  its 
p-norm,  and  |x|  :=  |x|2  is  its  Euclidean  norm. 


2.  BOUNDS  ON  THE  APPROXIMATION  ORDER 


2.1.  Principal  shift-invariant  spaces 

We  are  interested  in  characterizing  the  approximation  order  of  the  spline  spaces  {$h}h- 
This  is,  by  definition,  the  maximal  nonnegative  k  for  which 

distoo (f,Sh)  =  0(hk),  when  h  -*•  0, 


for  every  fc-admissible  /. 

In  order  for  this  definition  to  make  any  sense,  we  need  to  define  precisely  the  spaces  {s*}/u  as 
well  as  explain  the  notion  of  “^-admissible”.  We  start  with  the  former. 
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We  take  $h  to  be  an  appropriate  closure  of  the  linear  hull  of  the  /i2Zd-translates  of  some  ‘basic’ 
function.  Specifically,  we  take  s h  to  consist  of  functions  of  the  form  V’/iO  -  a)a(°0»  with, 

possibly,  some  restriction  imposed  on  the  coefficient  sequence  a.  Because  of  the  nature  of  the  results 
in  this  paper,  it  is  convenient  to  scale  up  Sh,  i.e.,  to  look  at  the  space 

(2.1)  &:={/(*•):  /€**}. 

The  space  Sh  is  a  principal  shift- invariant  space,  which  means,  by  definition,  that  it  is  “spanned” 
by  the  integer  translates  of  one  basic  function  fa  (which  happens  to  be  4>k(h‘)).  Denoting  by 

S(4) 

the  principal  shift-invariant  space  generated  by  the  integer  translates  of  <{>,  we  can  then  write 
Sh  -  S(fa).  Since  the  oo-norm  is  scale-invariant,  we  have 

dist  oo(/,s/i)  =  dist  oo(f(h-)>S(fa)), 

hence  the  change  from  Sh  to  the  scaled  space  Sh  =  S(fa)  requires  nothing  more  than  switching 
from  /  to  the  correspondingly  scaled  f(h-).  As  a  simple  example,  note  that  in  the  scaling  case, 
when  rfrh  is  the  h-dilate  of  fa ,  the  scale-up  procedure  undoes  the  dilation  and  hence  fa  =  fa  =  fa 
for  all  h.  In  other  words,  Sh  does  not  change  with  h. 

Thus  our  setting  is  as  follows:  we  hold  in  hand  a  collection  {Sh}h  of  spaces,  each  of  which  is  a 
principal  shift-invariant  space  generated  by  some  h-dependent  function  fa.  Then  for  a  “reasonable” 
function  /,  we  consider  the  quantities  dist00(/(h-),5/l).  Whenever  these  quantities  decay  to  0  like 
hk,  we  say  that  [Sh}h  provides  approximation  order  k  for  /.  If  dist** (f(h'),Sh)  =  0(hk)  for  all 
^-admissible  functions,  then  we  say  that  {Sh]h  provides  approximation  order  k. 

We  have  not  yet  defined  the  topology  used  in  the  definition  of  the  principal  shift-invdriant 
space  S(<f>).  While  the  derivation  of  lower  bounds  is  largely  independent  of  the  topology  used  in  the 
definition  of  this  “spline”  space  (since  only  a  small  subset  of  the  space  is  usually  employed  in  the 
analysis),  upper  bounds  are  intimately  related  to  the  way  S(<p)  is  defined:  results  on  upper  bounds 
become  stronger  with  the  weakening  of  the  topology  in  which  the  limits  )T)a€2Z*  <f>(-  -  o)a(a)  are 
calculated.  In  the  absence  of  a  standard  definition  for  the  space  S(cf>),  we  have  chosen  here  the 
following  one,  which  is  motivated  by  the  particular  way  in  which  we  shall  derive  upper  bounds  in 
the  next  section. 

Definition.  The  principal  shift-invariant  space  S(<f>)  is  the  space  of  all  locally  bounded  func¬ 
tions  fa' a,  for  which  the  double  sum 

fa\fa'd)  —  Yj  -  0)<KP  -  a)a(a) 
is  absolutely  convergent  for  every  x  €  Dl*1. 

If  <f>  has  compact  support,  then  S(<f>)  contains  fa' a  for  arbitrary  a.  Furthermore,  if  <f>  has  some 
decay  at  oo,  then  S(<f> )  contains  all  fa' a  for  which  a  does  not  grow  too  fast  at  oo.  Here  is  a  sample 
proposition: 
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Proposition  2.2.  Assume  that,  for  every  p  G  IIn,  the  series  converges  pointwise  absolutely 
to  a  locally  bounded  function,  and  let  An  be  the  space  of  all  sequences  a  :  7Ld  C  of  (at 
most )  polynomial  growth  n  at  oo.  Then  4>+'An  C  S(<f>).  In  particular,  (f>*'An  C  S(<j>)  in  case 
|^(x)|  =  0(|x|“m)  for  some  m  >  n  +  d,  as  x  -*  oo. 

Proof:  We  will  show  that,  for  a  G  An,  (<f>*'a)\xd  €  An,  from  which  it  will  follow  (because  of  the 
assumption  on  <f>)  that  converges  absolutely  to  a  locally  bounded  function,  and  therefore, 

by  the  definition  of  S(<f>),  <f>*'a  G  S(<f>).  Without  loss,  we  may  assume  that  both  <j>  and  a  are 
nonnegative  (otherwise  take  absolute  values). 

By  assumption,  we  can  find  a  constant  const  such  that  ||<£*'p||l00(c)  <  const  for  all  normalized 
monomials  xa/a!,  a  G  2Z+ ,  |a|i  <  n.  It  follows  that 

M+'pWl^C)  <  const  max  |ITp(0)| 
fr|i<» 

for  all  p  €  IIn.  Now,  let  y  G  IR1*,  and  set  y  =:  ty  +  ay,  with  ty  G  C  and  av  G  Z2d.  Since 
(<f>*'p)(-  +  oty)  =  <t>*'(p(-  +  Qfy)),  we  deduce  that  (<A*'p)(y)  is  the  value  at  ty  of  ^*'(p(.  +  ay)),  and 
therefore,  by  the  argument  above,  \<f>*'p(y)\  <  const  max|7|,<„  \D'1p(ay)\.  Thus  <f>*'p  =  0(|  •  |degp) 
at  oo,  and  hence  (<f>*'p)\xd  €  An  for  any  p  G  IIn. 

As  for  by  definition  of  An,  a  G  An  can  be  bounded  by  some  p  G  IIn  (in  the  sense  that 
a(a)  <  p(«)  for  all  a),  hence  is  dominated  by  <f>*'p  and  therefore  \xd  G  An.  4 

For  the  approach  taken  in  this  paper,  it  is  important  that  the  sum  <j>*'eg  be  well-defined  for  any 
exponential  e$,  0  G  HI*1.  Therefore,  we  assume  that  each  operator  <{>*'  is  well-defined  and  bounded 
as  a  map  from  l to  L^,  and  denote  the  corresponding  norm  by  Some  conditions  related 

to  the  boundedness  of  ||<£*'||  are  recorded  in  the  following  proposition. 

Proposition  2.3.  The  norm  of  the  operator  <p*'  is  ||  Ylaen*  “  a)l  ll>  hence,  this  operator  is 
bounded  if  and  only  if  the  series  Yja&L*  | <(>(•  -  a)|  is  pointwise  convergent  to  a  bounded  function. 

Proof:  For  any  a  G  1 00  and  any  x  G  Dld, 

(2.4)  |(^*,o)(x)|  <  Wx  “  *)l  NL<*->» 

ae.HL* 


which  shows  that  \\<f>*'\\  <  ||  Y.a&L*  “  Q)l  II- 

For  the  converse,  let  K  be  a  finite  subset  of  2Zd  and  let  E  be  the  collection  of  all  sequences 
e  :  7Ld  -*  {0,±1}  which  are  supported  on  K.  For  every  x  G  IRd,  there  exists  e  G  E  such  that 
=  £o6k  \4>{x  -  o)|.  Since  E  is  finite,  we  conclude  that 

M  >  *»p  Ill'll  >  II Y,  «•  -  «)l  II- 

*€B  a£K 

Writing  ZZd  as  the  union  of  a  nested  collection  {A'}  of  finite  sets,  and  using  the  fact  that 

11 E  w-  -  °)i  11  — .  11 E  w-  -  «)i  11. 

a£K  K~*1L  aen* 
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we  obtain  the  desired  result.  4 

The  last  proposition  implies  that  (f>  £  Li(TRd)  whenever  <f>*'  is  bounded,  and,  hence,  that  the 
Fourier  transform  <j>  of  <j>  is  a  well-defined  continuous  function.  Also,  a  sufficient  condition  for  the 
boundedness  of  <f>*'  is  the  integrability  of  the  maximal  function  <£#(x)  :=  ||^||l00(i+c)- 


2.2.  Admissibility 

Next,  we  turn  to  the  definition  of  the  space  of  admissible  functions  associated  with  the  oo-norm: 

Definition.  A  function  f  of  at  most  polynomial  growth  at  oo  is  termed  here  ^-admissible  if 
(1  +  |  •  |fc)/  is  a  Radon  measure  of  finite  total  mass.  For  such  a  function  /,  we  denote  by 

11/111 


the  total  mass  of  (1  +  |  •  | *)/. 

It  follows  that  /  is  fc-admissible  (for  some  k  >  0)  only  if  /  is  a  measure  of  finite  total  mass. 
In  particular,  any  admissible  /  is  bounded.  It  is  worthwhile  to  keep  in  mind  two  examples  of 

j  A 

admisssible  functions.  The  first  is  the  exponential  f  =  e$,  9  £  JR  .  In  this  case,  /  =  and  since 
/  is  compactly  supported,  /  is  admissible  of  all  orders.  However,  ||/||^  =  1  +  \9\k,  and  this  grows 
with  k  and/or  9.  The  other  example  occurs  when  /  is  a  function.  In  this  case,  /  is  A:-admissible 
whenever  (1  -f  |  •  | k)f  £  L\{ IRrf). 

As  usual,  in  case  k  is  integral,  the  admissibility  condition  can  be  interpreted  in  terms  of  the 
fcth  order  derivatives  of  /: 

Proposition  2.5.  A  function  f  is  k-admissible  for  some  k  £  ZZ+  if  and  only  if  the  Fourier  trans¬ 
forms  of  f  and  of  all  its  kth  order  derivatives  are  measures  of  finite  total  mass. 

*****  w 

Proof:  Let  fa  be  the  ath  order  (distributional)  derivative  of  /,  hence  fa  :w  i-+  (iti>)“/(w),  and 

choose  ca  so  that  £oc<»ll0,l  =  M*  for  x  €  IRrf,  i.e.,  cQ  =  (£)  for  |a|i  =  k  and  ca  =  0  otherwise. 

Then  ]£acal/«l  =  j  •  |*!/|,  therefore,  if  fa  is  a  measure  of  finite  total  mass  for  each  a  £  ZZ+  with 

|a|j  =  fc,  then  so  is  |  •  |*|/|,  hence  so  is  |  - 1* \f\.  Thus,  if  also  /  is  of  finite  mass,  then  we  conclude 
,  ^  *** 
that  so  is  (1  +  |  ■  |  )|/|.  The  converse  is  even  simpler:  if  /  is  ^-admissible,  then  /,  as  well  as 

A  ,  A 

w  h-»  ( iw)°f(w )  for  |a|i  =  k,  are  majorized  by  a  measure  of  finite  mass  (viz.  (1  + 1  •  |k)/),  and  hence 
the  Fourier  transform  of  /  and  of  all  its  derivatives  of  order  k  are  measures  of  finite  mass.  A 


2.3.  Upper  bounds 

We  obtain  upper  bounds  for  the  approximation  order  by  considering  approximation  to  expo¬ 
nentials  e$,  6  €  HI*1.  Our  starting  point  is  the  following  result  from  [R2]: 

Result  2.0.  Let  9  £  VHd,  and  assume  that  the  sequence  {<f>h}h  satisfies  the  following  conditions: 

(a)  suppfa  C  B,  for  all  h,  and  for  some  h-indcpendcnt  compact  B. 

(b)  The  functions  {<f>h)h  are  uniformly  bounded. 
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Then, 

(2.7) 


\\<f>h*'ehe  -  4>h{h6)eh9\\  <  cdist00(e/,tf,5(^)). 

The  proof  provided  here  will  make  use  of  the  following  condition,  which  is  a  consequence  of 
(a)+(b),  but  implies  only  (b): 

(ab)  sup/j  ||<^/i*,||  <  oo. 

Proof:  Fix  h,  and  let  /  €  S(<j>k)-  Since  <f>h*'g  =  g*'<f>h  for  all  g  €  S((j>h),  by  [Bl],  and  also 
e9*'4>h  =  for  any  d,  we  have 

(2.8)  || <t>h*'eM  -  4>h{hO)eM\\  <  \\<J>h*'eh8  ~  4>h*' f\\  +  \\f*'<i>K  -  ehe*'4>h\\  <  2||<^*'||  || eh9  -  /||. 

Since  /  €  S(<j>h)  was  arbitrary,  the  result  follows,  with  c  =  2supA  ||^*'||.  4 

Since  the  key  to  the  above  argument  is  the  “flip”  property:  <f>+' f  =  €  S(<f>),  we 

can  extend  the  result  to  any  <f> ,  compactly  supported  or  not,  with  that  property.  Our  particular 
definition  of  the  space  S(<f>)  was  chosen  primarily  to  ensure  the  “ Hip ”  property. 

Flip  Lemma  2.9.  For  every  f  6  S(<f>), 


Proof:  The  argument  follows  the  one  given  in  [Bl].  We  fix  /  €  S(<f>)  and  z  €  IRd,  and  wish  to 
show  that  both  <f)*'f(x)  and  f*'<f>{x)  converge,  and  to  the  same  limit.  Assuming  /  =  <f>*'a  for  some 
sequence  a,  we  write  explicitly 

(^*'/)(x)  =  <f>(x  -  a)  £  4>{a  -  /?)a(/3). 

ae.iL*  pen* 

By  the  definition  of  S(<f> ),  this  double  sum  is  absolutely  convergent,  hence  we  may  rearrange  terms 
(and  replace  a  by  a  -f  /?)  to  get: 

E  #«)  E  «*-(«+«w) =/*'#*)• 

aen*  pen* 


* 

Theorem  2.10.  Assume  that  the  operator  <f>*’  is  bounded.  Then,  for  any  0  6  md, 

(2.11)  \\<t>*'eM  -  t(hO)eM ||  <  2||**'||  dist oo(ew,  $(*)). 


Proof:  Repeat  the  proof  of  Result  2.6,  but  replace  the  reference  to  [Bl]  by  a  reference  to  Lemma 
2.9.  6 

Theorem  2.10  can  be  interpreted  in  two  different  ways.  On  the  one  hand,  it  suggests  that  a 
‘near-optimal’  approximant  for  the  exponential  e^t  from 

Sh  :=  S(*a) 

is  provided  by  (the  supposedly  well-defined)  <f>h{hO)~x<t>h*'eM-  The  following  corollary  records  this 
fact. 
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Corollary  2.12.  Assuming  that  fa*1  is  bounded  and  that  fa{hO )  ^  0,  we  get 


(2.13) 


\\fa(h0)~lfa*'ek6 


-«m||  <2 


\\M 

\Uh9)\ 


distoo(e/i0, 5/i). 


We  note  that  the  ratio  \\fa*'\\/\<j>h(h0)\  is  independent  of  the  way  cj>h  is  normalized,  and  hence, 
the  right  hand  side  of  (2.13)  is  independent  of  the  particular  normalization  we  choose  for  <f>h-  But 
the  estimate  (2.13)  is  useful  for  the  derivation  of  bounds  for  the  approximation  order  only  in  case 
the  sequence  {\\<t>h*'\\/<t>(h0)}k  is  bounded. 

Fortunately,  Theorem  2.10  can  be  used  directly  to  derive  upper  bounds.  We  simply  observe 
that,  in  case  the  operators  {<f>h*'}h  are  uniformly  bounded,  Theorem  2.10  shows  that  distoo(e/,«,5/i)  = 
0(hk)  only  if  the  same  holds  for  || fa*'eh$  ~  ^h(h0)eile\\,  i.e.,  for  \\e^ha(fa*'chs)  -  fa(h0)\\  (since 
|cas(x)|  =  1  for  every  9  €  IR**  and  every  x  £  IR**).  Since  e-heifa^Pht)  =  (e_^^)*'l,  we  obtain 

(2.14)  ||(c-/»«^fc)*,l  -  MM)\\  <  cdist00(efc«,<S!/i). 


Thus,  if  we  assume  that  we  have  approximation  order  k,  then  we  must  have 

(2.15)  ||(«-M*fcXl  -  fa(M) ||  =  0(hk). 

Since  the  function  (e^hafa)*1 1  is  2Zd-periodic,  (2.15)  implies  that  its  Fourier  coefficients  (excluding 
the  O’th  coefficient)  must  be  of  size  0(hk).  Furthermore,  (2.15)  implies,  in  particular,  that 

(2.16)  IRe-w^K1  ~  ?fc(M)IUa  (c)  =  0{hk), 

which  means  that  the  2-norm  of  the  Fourier  coefficient  sequence  for  this  periodic  function  is  of 
order  0(hk).  Since  fa{h9)  is  part  of  the  constant  term  of  this  function,  these  coefficients  are 
(«-M**)*'ir09)  for  P  €  2x7Ld\0.  We  compute 

({e-hethb'lTiP)  =  f  JZ  e~M “  a)^(*  “  a)e-p( 0  dt 

(2.17)  =  E  /  <-w(0^(0«-/»(0 dt 

a&t*  JC~a 

=  /  fa{t)C-hS-p(t)  dt  =  +  P)i 

J  R' 

where,  for  the  second  equality,  the  fact  that  fa  £  L\  was  used.  Therefore,  we  conclude  that 

As  a  matter  of  fact,  nothing  in  the  above  arguments  requires  the  approximation  order  to  behave 
like  a  power  of  h ,  and  we  thus  obtain  the  following. 
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The  Upper  Bound  Theorem  2.18.  Assume  that  the  fa*'  are  bounded,  and  let  0  €  lRrf.  If 

dist  (cm,  5(^/0)  =  0(pe(h)) 
for  some  (univariate)  function  p$,  then ,  for  every  h, 

Yl  I  MM  +  0)\2  <  constH^t'Hp^/i)2. 

pe.2itHd\o 

In  particular,  if  we  normalize  {fa}h  to  obtain  a  uniformly  bounded  [<f>h*'}h,  then  {S(<j>h)}h  provides 
approximation  order  k  to  the  exponential  function  eg  ,  6  6  IRd,  only  if 

£  \Mhe+p)\2<c6h2k. 

P&*TZ*\0 

Note  that  the  above  implies  that,  in  order  to  obtain  ^-approximation  order,  it  is  necessary  to 

have 

(2.19)  | fa(h0  +  P)\  <  c8hkt  0  €  2*^0,  0  €  IR*'. 

It  is  this  slightly  weaker  condition  that  we  use  throughout  the  paper  in  order  to  obtain  upper 
bounds  on  the  approximation  orders. 

It  is  remarkable  that  the  result  avoids  an  application  of  Poisson's  summation  formula  (namely, 
the  convergence  of  the  Fourier  series  of  (e_A®0/,)*'l  was  not  required),  and  hence  no  smoothness 
conditions  were  imposed  on  {<f>h}h-  Also,  no  “regularity”  condition  was  needed  in  the  upper  bound 
theorem,  i.e.,  neither  {fa(ty}h  nor  {^(0)}A  were  required  to  stay  away  from  0.  (However,  it  is 
plausible  that,  in  the  singular  cases,  this  upper  bound  overestimates  the  actual  approximation  order 
by  the  order  of  the  zero  of  h*-*  <f>h(h0)  at  h  =  0.) 

The  upper  bounds  were  derived  under  the  assumption  that  the  exponential  function  e#  is 
admissible,  hence  should  be  approximated  well.  Under  a  stronger  assumption  on  the  rates  of  decay 
of  each  fa  at  oo,  we  can  show  that  the  Upper  Bound  Theorem  2.18  remains  valid  even  if  we  only 
wish  to  approximate  the  test  functions  in  V,  i.e.,  infinitely  smooth  compactly  supported  functions. 

Theorem  2.20.  Assume  that  {||^/i*'||}*  is  bounded  and,  in  addition,  {fajh  satisfies  the  condition 

(2.21)  £  \fa(x  -  o)|  <  chk,  xeC. 

\o\>l/h 

Then  {S(fa)}h  provides  approximation  order  k  to  all  functions  in  V  only  if,  for  every  0  €  Hf*  and 
every  0  €  2xZZ<<\0, 

\fa(hf  +  fi)\  <c,hk, 

where  e#  is  independent  of  0. 

Proof:  We  fix  9  €  Dlrf,  and  choose  /  €  V  with  ||/||  =  1  such  that  /  =  e#  on  some  neighborhood 
of  0,  e.g.,  on  3C.  Let  x  €  C,  h  <  1,  and  set  A  :=  /(/»•).  Then,  since  |o|  <  1  fh  implies  that 
C-aC  3  C/h,  we  have  fh(  x  -  o)  =  e*#(z  -  o)  for  |o|  <  l/h,  and  therefore 

l(A -<«)•>»<*)!  <s  £  IM«)l  =  0(*‘). 

\o\>l/h 
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In  a  similar  fashion,  | <f>h*'Uh  -  cm)(x)|  =  0(hk),  too. 

As  in  Result  2.6,  we  obtain  that,  for  each  g  6  S(<j>h ), 

Mh*'A  -  h*% ||  <  || -  5)11  +  I ||(A  -  9)*'+k II  <  2||<A*'||  II A  -  g\\. 

Therefore,  for  x  €  C  we  obtain  that 

\(e-he<f>h)*'l(x)  -  &,(/i0)|  =  | <t>h*'ehe(x)  -  ^a(M)«A«(*)|  =  I <fo*'e/»«(x)  -  «a«*'^(x)| 

<  |4**'A(*)  ~  A*'^(*)l  +  0(hk)  <  2||^*'||  distooCA,  Sh)  +  0(hk). 

Since  the  function  c_as<^a*,1  is  Z^-periodic,  and  since  we  assume  the  boundedness  of  {<f>h.*'}iu  we 
conclude  that 

-  4>{h0) ||  <  const distooCA,  5a)  +  0(hk). 

Finally,  in  the  prov.f  of  the  Upper  Bound  Theorem  2.18  we  have  observed  that  {^(/ifl-f /S)}^,^  are 
the  Fourier  co'  f’kients  of  t’.  Unction  (c-a«^a)*,l»  and  therefore,  if  {5a}a  provides  approximation 
order  k  for  /,  we  vrst 11  we 

\4(he  +  p)\  <chk,  /?€2 ttZZ^O, 

where  c  depei  :  -a  /,  i.e.,  on  0,  but  is  independent  of  (3 .  4 


Lower  bounds 

The  ’ippir  bound  anal' sis  provides  us  with  near-optimal  approximants  £h$  from  5a  to  the 
xponential  g.  tyt  u  closely  related  approximants  to  establish  lower  bounds  on  the  approxima¬ 
tion  orders  province  oy  {5a}a<  The  idea  of  producing  an  approximant  from  5a  for  the  function 
A  :=  f(h-)  is  very  si  aple:  we  write 

A(*)  =  Wd  f  f(0)eh$(x)  d9 , 

Jr4 

hence  can  provide  the  approximant  to  A  in  the  form 

(2jr)~rf  f  f(9)ehe(x)  d9. 

Jn4 

Of  course,  we  still  need  to  make  sure  that  this  approximant  lies  in  5a,  and  that  /  is  “reasonable” 
enough  for  the  above  integration  to  make  sense.  But,  even  in  such  a  case,  the  above  approximation 
scheme  is  “too  global”  (in  the  sense  that  all  the  Fourier  transform  information  of  /  is  taken  into 
account),  and  therefore,  in  order  to  simplify  our  error  analysis,  we  use  a  suitable  nonnegative 
continuous  cut-off  function  a  with  support  near  the  origin.  For  convenience,  we  assume  that  supper 
lies  in  the  ball  Bv  -  {x  €  IRd  :  |x|  <  »?},  that  a  is  1  on  Bv/2,  and  that  ||er||  =  1.  Furthermore, 
since  the  approximation  scheme  should  be  applicable  to  any  admissible  function  /,  we  only  know 

A  A  A 

that  /  is  a  measure,  and  therefore  prefer  the  notation  df(9)  to  the  notation  f(9)  d9.  Thus,  our 
approximation  scheme  takes  the  form: 

(2.22)  A(x)  :=  fihx)  w  Rhf{x)  :=  (27r )— /  ehe(x)cr{h9)df{9). 

Jn 4 
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Then 


(2.23) 


(2x)*\Rhf-fh\<  (  \eM<T(h9)-ehe\\df\(6) 

Jr- 


<  f  KM)  -  1|  |i/|(*)  +  /  ff(M)|£M/e„  -  1|  |^(«). 

Jnd  Jnd 

If  now  /  is  ^-admissible,  then 

(2A/v)-*  fa-  «(m»  i  am  <  f  w*(i  -  cm)  i<tfi(»)  rt 0 

JWi*  JRd  h~'° 

(since  the  integrand  vanishes  on  5„/(2A)),  and  therefore  the  first  integral  in  the  last  line  of  (2.23) 
is  o(hk).  For  the  second  integral,  we  need  to  look  more  carefully  at  the  ‘periodized’  error 


(2.24) 


Chs/ehe  - 1. 


2.4.  Lower  bounds:  analysis 

We  choose  the  ‘near  optimal’  approximation  e he  from  Sh  to  e^i  in  the  form 


(2.25) 


„  <t>h*'ehe 
£ho  •“  ^ 


MhO) 

With  this  choice,  the  approximation  Rhf  takes  the  more  explicit  form 
(2.26)  Rhf  := 


with  fH  the  bounded  analytic  function 

(2.27)  /{(*)  *  (?n)~d  f  4 

J  R*  Oi 


ihx-9 


-o{h0)  df{6). 


<f>h{hO) 

In  particular,  this  makes  clear  that  Rhf  is  indeed  an  element  of  Sh- 
Further,  the  ‘periodized’  error  takes  the  form 


ulu  $hm 


Recall  from  (2.17)  that 


((e-fc^fc)*'ir(/J)  =  MhO  +  /?),  fie  2  *7Zd. 
Consequently  the  Fourier  series  of  the  periodized  error  has  the  form 

y'  <i>h(hO  +  P) c  _  ^  _  V'  ^h(h0  +  fi) ^ 

/3€2ir  TL*  4>h{h8)  /3€27rZS<*\0  ?h{hO) 


and  it  always  converges  to  the  periodized  error,  at  least  in  a  distributional  sense.  While  the 
estimates  for  this  sum  provided  by  summability  methods  seem  to  be  hard  to  analyze,  we  can 
expect  that,  for  a  smooth  <f>h,  the  series  will  converge  absolutely.  In  such  a  case,  we  obtain  the 
following  important  estimate. 
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Theorem  2.28.  For  an y  h,  for  which  fa*'  is  bounded, 


(2.29) 


I \sh§/ekt  -  1||  =  || 


k(hB) 


e/ioll  < 


\$h(M  +  0)\ 
0€2»:Z<*\O  i^(M)l 


It  may  seem  surprising  that  we  did  not  use  here  the  approximation 


(2.30) 


CM  « 


<j>h*'eh9 

km 


derived  during  the  discussion  of  upper  bounds.  The  reason  is  simple.  Recall  from  (1.1)  that 
$h(h8)  =  (<M'cm)(0)  =  (c_m^/i*'1)(0),  and  this  equals  T,pe2*7td(e-M<l>h*'lT(P)-  Therefore,  by 
(2.17),  $h(h9)  =  $h(h9)  +  E^62JrZ-'\o  +  £)•  Hence  one  of  the  two  sums 

V'  -h  P)\  an(j  y'  I  $h{h9  +  /?)[ 

0€2»Z'\O  /362irZ^\0 


goes  to  zero  with  h ,  then  they  go  to  zero  at  exactly  the  same  rate,  since  lim^—o  <j>h{M)/<i>h(h9)  =  1 
in  such  a  case.  This  means  that  we  lose  nothing  in  the  estimate  (2.29)  if  we  use  the  approximation 
ehi  =  instead  of  but  gain  simplicity,  since  the  Fourier  transform  4>h  is  usually  more 

4k(ht)  ^*(M) 

readily  accessible  than  the  symbol  4>h. 

In  the  sequel  we  exclusively  use  the  right  hand  side  of  (2.29)  to  bound  ||£m/cm  -  1||>  hence 
obtain  positive  approximation  orders  only  when  the  right  hand  side  of  (2.29)  tends  to  0  with  h. 


2.5.  Lower  bounds:  synthesis 
With  the  bound  (2.29)  in  hand,  we  return  to  the  error  estimate 

(2.31)  (ix)d\Rhf-fk\<  I  (1  -  a(M»  |d7l(0)  +  /  *(M)I«»./«m  -  11  \i?\W\ 

J  iv4  J  nd 

see  (2.23).  Having  observed  earlier  that  the  first  integral  is  o(hk)  whenever  /  is  ^-admissible,  we 
now  want  to  make  the  second  integral  0(hk).  By  Theorem  2.28,  we  have  in  hand  the  estimate 


(2.32)  /  a(h8)\eM/eM  -  1|  \df\(8)  <  [ 

J  JB 


T 

Bi/k  /J62jrErf\0  l^»(^)l 


\Mh8  +  f3)\ 


We  know  from  the  upper  bound  discussion  that  approximation  order  k  requires  that 
0(hk)  for  fixed  8  uniformly  in  /?.  This  suggests  the  assumption  that 


I _ 

|*fc(M)l 


ii  E 

P&*TLd\0 


i  Uh-m 
\Uh-)\ 


Ikoote,)  =  0(hk). 
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Unfortunately,  except  for  the  case  of  spectral  approximation  orders  (see  Section  3.3),  not  many 
examples  of  interest  satisfy  this  assumption.  We  choose  instead  the  following  more  subtle  condition, 

in  which  we  employ  functions  Vh(x)  to  describe  the  behavior  of  ^$1-— ' for  i  near  0.  In  fact,  for 

IM*)I 

any  indexed  set  of  positive  functions  on  2?„,  we  have 


(2.33) 


L 


£ 

B*Jk  /J€2*®*\0 


\MM  +  P)\ 
\Uh0)\ 


me)-  £  %.<*,>  L/l(M)l<,/1w- 

P£2ir7li\0 


Consequently, 


(2.34) 


/  o(h0) \ehe/ehe  - 1|  \df\(B)  <  A(u,  tj)  f 
JRd  Jb 


BJh 


vh(hQ)  |d/|(0), 


with 


T  £ 

/?  £2  ttZS*4  \0 


httei  I. 


The  Synthesis  Condition.  We  say  that  {4>h}h  satisfies  the  {/-synthesis  condition  if,  for  some 
indexed  set  {{//,}/,  of  functions  defined  in  some  neighborhood  of  0,  and  some  rj  >  0,  A(v,ri)  is  finite. 
This  means  that  the  number 


(2.35)  A(v)  :=  inf  A(u,  rj) 

v 

is  finite. 

Now  notice  that  /R,<(|h0|fc  +  hk )  |d/|(0)  =  A*||/||'fc.  This  suggests  the  choice 

(2.36)  i/fc(*)  :=  \x\k  +  hk. 

We  call  the  {/-synthesis  condition  with  respect  to  this  u  the  synthesis  condition  of  order  k. 
With  this,  we  infer  from  (2.31)  and  (2.34)  the  following: 

The  Lower  Bound  Theorem  2.37.  Assume  that  each  4>h*'  is  bounded,  and  that  {<f>h}h  satisfies 
the  synthesis  condition  of  order  k.  Then,  for  every  k-admissible  f,  ||  Rhf  -  /(/i-)||  =  0(hk),  with 
Rhf  defined  as  in  (2.22).  In  particular,  [S(4>h.)}h  provides  approximation  order  k. 

As  a  matter  of  fact,  the  above  discussion  provides  alsu  the  following  significant  information  on 
the  constant  in  the  0(hk)  expression. 

Corollary  2.38.  Assuming  the  synthesis  condition  of  order  k,  we  have,  for  every  k-admissible  f, 

n^*7;  -  aii  <  M-w  h/imm + 

with  A(u)  as  defined  in  (2.35)  (for  Uh(x)  :=  |x|*  +  hk). 

Stronger  results  can  be  obtained  under  more  restrictive  assumptions  on  the  functions  u  in  the 
synthesis  condition.  For  example,  if  Uh{x)  is  independent  of  x  (but,  of  course,  depends  on  h ),  then 
(2.3!)  and  (2.34)  allow  us  to  conclude  the  following  improved  version  of  the  last  corollary: 
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Corollary  2.39.  If  {<j>h}h  satisfies  the  u- synthesis  condition ,  with  each  Vh  being  a  constant  func¬ 
tion,  then,  for  every  k-admissibie  function  f, 

«**•'/»  -  All  <  (2i)-V/.(0)||/M^)  +  «(**). 

In  particular,  if  t//,(0)  =  o(hk),  then  so  is  dist00(//j,5/l). 

3.  Applications 

We  apply  here  the  general  results  of  the  previous  section  to  the  three  main  families  of  basic 
functions:  basic  functions  obtained  by  differencing  another  basic  function,  h-independent  basic 
functions,  and  basic  functions  which  are  compactly  supported  uniformly  in  h.  These  families 
overlap,  and  many  specific  examples  fall  into  two  or  three  of  these  categories  (e.g.,  polynomial  box 
splines  satisfy  all  of  the  above  conditions). 


3.1.  The  scaling  case 

We  •  H1  th'1  terminology  “the  scaling  case”  if  the  functions  {<f>h}h  are  actually  independent  of  h. 
This  m«  *  ,  •  that  all  the  spaces  Sh  =  S(<f>h)  are  the  same,  or,  in  other  words,  that  the  original  spaces 
{sh}h  are  obtained  by  dilating  S(<f>)  (with  <j>  :=  <f>i,  say).  In  this  case,  it  is  convenient  to  speak  of 
the  approximation  order  provided  by  <f> •  and  mean  by  this  the  approximation  order  provided  by 
{Sh.  -  S{(f>)}h •  The  upper  bound  theorem  specializes  in  this  case  to  the  following. 

Theorem  3.1.  Let  <f>  be  a  measurable  function  whose  associated  operator  4>*'  is  bounded,  and 
whose  Fourier  transform  is  k  times  differentiable  at  every  f3  €  2n7Ld\0  for  some  k  €  2Z+.  Then  <f> 
provides  approximation  order  k  only  if  <f>  has  a  zero  of  order  k  at  every  /3  €  2ir2Z  \Q. 

Proof:  Since  <t>h~  4>  for  every  h,  we  must  have,  by  the  Upper  Bound  Theorem  2.18, 

${h6  + 13)  =  0{hk),  V0eIRd. 

A 

By  taking  sufficiently  (but  finitely)  many  0’s,  we  conclude  that  <f>  must  have  a  zero  of  order  k  at 
each  /?  €  27rZd\0.  4 

The  scaling  case  was  analyzed  in  great  detail  in  the  literature,  especially  for  a  compactly 
supported  <f) ,  and  the  various  results  for  this  case  are  usually  in  terms  of  “Strang- Fix  Conditions”, 
[SF].  We  record  below  the  following  version  of  these  conditions,  whose  sufficiency  follows  from 
[DM1]  and  [Bl],  while  their  necessity  was  proved  in  [R2]  (and  previously,  in  a  weaker  form,  in 
[BJ]). 

Result  3.2.  Assume  that  <j>  is  a  bounded  measurable  compactly  supported  function  that  satisfies 
1(0)  0,  and  let  k  €  IN.  Then,  4>  provides  approximation  order  k  for  every  f  in  the  Sobolev 

space  W £,  if  and  only  if  II<fe  C  S(<f>). 

Exploiting  the  Upper  Bound  Theorem  2.18,  we  are  able  to  extend  the  “only  if’  statement  in 
Result  3.2  as  follows: 
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Theorem  3.3.  Let  0  be  a  bounded  mensurable  compactly  supported  function.  Then  the  ap¬ 
proximation  order  provided  by  <f>  to  functions  in  V  cannot  exceed  the  degree  of  the  least  degree 
polynomial  p  that  satisfies  <j>*'p  £  II. 

Proof:  By  Theorem  2.20,  the  Upper  Bound  Theorem  2.18  holds  here  (although  only  approxima- 

A.  A 

tion  to  functions  in  V  is  required).  Thus,  by  Theorem  3.1  (and  the  fact  that  0  is  entire),  0  must 
have  a  zero  of  order  k  at  each  /3  6  2xEd\0,  in  case  the  approximation  order  is  k.  Now,  it  is  known 
(cf.  e.g.,  [BR])  that,  for  a  compactly  supported  distribution  0,  0*' II<jt  C  II  if  and  only  if  <f>  has  a 
zero  of  order  k  at  each  (3  €  27rZZd\0.  Thus  the  degree  of  the  least  degree  polynomial  p  for  which 
<f>*'p  #  II  is  an  upper  bound  for  the  appoximation  order.  4 

In  case  0(0)  ^  0,  the  above  result  reproduces  the  “only  if’  statement  of  Result  3.2,  since  it  is 
well-known  (cf.  [Bl])  that  in  such  a  case  <f>*'p  6  II  if  and  only  if  p  6  5(0). 

The  above  result  does  not  admit  a  direct  extension  to  the  case  of  global  support,  since  in 
general  the  operator  0*'  need  not  to  be  defined  on  polynomials  in  that  case.  However,  we  always 
have  the  following: 

Theorem  3.4.  Assume  that  0 *'  is  bounded.  Then  the  approximation  order  provided  by  0  is 
positive  only  if  0+'l  =  const. 

^  , 

Proof:  By  Theorem  3.1,  0  vanishes  on  2?r2Z  \0.  Let  {j)n}  be  an  approximate  identity,  and  define 
rn  :=  r}n  *  0.  Then,  for  every  n,  t„  vanishes  on  2ir7Ld\0 ,  and  hence,  by  Poisson’s  summation 
formula, 

r„.'l=  £  ?.(/>)«„  =  ?(0). 

/J€2*ZJ 

A  A 

Thus,  r]n  *  (0*1)  =  0(0)  for  every  n,  and  by  letting  n  -*■  oo,  we  obtain  0*'l  =  0(0).  6 

Similar  results  can  be  obtained  with  respect  to  higher  degree  polynomials,  K  we  assume  that 
0*'  is  well-defined  on  such  polynomials.  We  omit  these  details  here. 

A. 

We  now  turn  to  the  synthesis  condition,  and  examine  the  quotient  In  many  of  the 

examples,  this  ratio  can  be  factored  into  two  terms 

=  *(*  + 

0(x) 

where  K(x )  decays  fast  enough  at  oo,  while  H(x)  has  a  zero  at  the  origin  of  high  enough  order.  In 
such  a  situation,  we  obtain  the  following  result. 

A 

Theorem  3.5.  Suppose  that  0*'  is  bounded  and  0(0)  ^  0.  Suppose  further  that,  for  some 
continuous  K  and  smooth  H, 

S|±^l  =  A'(x  +  /3Wx), 

0(l) 

and  that,  for  some  tj>  0, 

p£UKi\  0 
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If  K((3)  ^  0  for  some  (3  €  2ir7Ld\Q,  then  the  approximation  order  provided  by  <j>  equals  the  order 
of  the  zero  of  H  at  the  origin. 

Proof:  Let  (3  6  2it7Ld\0  be  such  that  K{(3)  0.  Let  k  be  the  exact  order  of  the  zero  of  H  at  the 

origin.  Since  H  is  smooth  at  0,  there  exists  9  €  IRd\0  such  that  H(h9)  ^  o(hfc),  and  consequently, 


<f>(h6  +  (3) 
$(h9) 


=  K(h9  +  (3)H(h6)  ^  o(hk). 


Since  tf>*'  is  bounded,  we  can  appeal  to  the  Upper  Bound  Theorem  2.18  to  conclude  that  the  relevant 
approximation  order  is  at  most  k. 

Next,  we  show  that  the  synthesis  condition  of  order  k  holds:  by  the  assumption  on  H ,  |  •  \~kH 
is  bounded  in  a  neighborhood  of  the  origin.  Thus,  for  a  small  enough  rj  and  any  16  5,, 


E 

0£2ir7Zd\Q 


\K(x  +  {3)H(x)\ 


peuTt^o 


This  shows  that  the  synthesis  condition  of  order  k  holds  here,  and  hence,  by  the  Lower  Bound 
Theorem  2.37,  the  approximation  order  is  at  least  k.  4 


3.2.  The  differencing  case 

Almost  all  the  basic  functions  now  in  the  literature  belong  to  this  family.  Here,  one  starts  with 
a  function  4>  which  has  polynomial  or  even  exponential  growth  at  oo,  but  for  which  p(D)<j>  decays 
at  oo  for  some  linear  differential  operator  p(D)  with  constant  coefficients  (e.g.,  <f>  is  a  fundamental 
solution  of  p(D)).  This  means  that,  away  from  the  zero  set  of  p(i-),  4>  coincides  with  K/p(i‘),  for 
some  smooth  function  K.  To  obtain  from  <t>  a  function  in  (a  suitable  closure  of)  S(<f>)  which  decays 
at  oo  (namely,  to  localize  <j>),  or,  equivalently,  to  remove  the  singularities  of  <f> ,  one  approximates 
the  differential  operator  p(D)  by  a  (finite-  or  infinite-)  difference  operator  T  which  is  supported 
on  2Zd,  i.e.,  approximates  the  polynomial  p(i*)  by  some  periodic  function  u.  In  order  to  make 
this  process  feasible  and  as  simple  as  possible,  the  real  variety  of  the  polynomial  p(i-)  should  be 
extremely  simple.  For  example,  in  the  box  spline  case  p  is  chosen  as  a  product  of  (more  or  less 
arbitrary)  linear  polynomials,  while  for  typical  radial  basis  (and  related)  functions,  the  operator 
p(D)  is  elliptic  (for  this  reason  we  will  refer  to  the  latter  class  of  basic  functions  as  belonging  to 
“the  elliptic  case”),  hence  one  has  only  to  resolve  the  singularity  of  <f>  at  the  origin.  For  h  >  0,  one 
replaces  p(D)  by  the  operator  ph{D)  :=  hde$pp(h~1D),  associates  Ph{D )  with  a  function  <j>h  in  a 
way  analogous  to  the  case  h  =  1  (so  that  now  Ph{i’)$h  =  K,  for  the  same  K  as  before),  and  repeats 
the  differencing  process.  In  case  p  is  homogeneous,  p/,  does  not  change  with  h,  hence  <ph  =  0,  all  h. 

For  box  splines,  one  easily  obtains  in  this  way  a  compactly  supported  function  T<£  ([BH],  [Rl]). 
In  contrast,  in  the  (homogeneous)  elliptic  case,  the  major  effort  was  devoted  to  the  localization 
process,  aiming  at  constructing  T  in  such  a  way  that  (T <f>y  be  as  smooth  as  possible  and  (T 1 
have  a  high  order  zero  at  the  origin.  The  simple  approximation  scheme  then  suggested  (cf.  e.g., 
[Jj,  [Bul-3],  [DJLR],  (P))  is  f(h-)  m  T <f>*’ f{h-),  and  the  convergence  rate  is  0{hk\\ogh\)  where 
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k  depends  on  the  rate  of  decay  of  T <f>  at  oo  and  on  the  order  of  the  zero  of  (T<py~  1  at  the 
origin,  but  never  exceeds  degp.  Under  further  assumptions  on  the  decay  rates  of  T<p  at  oo,  the 
approximation  rates  were  improved  to  0(hk),  but  in  some  cases  it  was  proved  that  0(/ifc|log/i|)  is 
the  exact  order  of  this  approximation  scheme.  None  of  the  above  references  provides  upper  bounds 
on  the  approximation  order.  Details  for  specific  cases  are  given  in  the  next  section. 

We  show  now  that  the  exact  approximation  order  in  all  elliptic  cases  in  the  literature  is  the 
degree  of  the  underlying  differential  operator  p(D),  regardless  of  the  decay  rates  ofT<p  at  oo,  or 
the  smoothness  of  (T <py~  at  the  origin. 

Theorem  3.6.  Let  p(D )  be  a  homogeneous  elliptic  operator  of  order  k  >  d,  d  being  the  spatial 
dimension.  Assume  that  <p  satisfies  the  equation 

P(  *-)$  =  K, 

where  K  is  a  continuous  bounded  function,  A'(0)  ^  0,  and  K(/3 )  ^  0  for  some  /?  €  27rZd\0. 
Assume  that,  for  some  sequence  c  of  polynomial  growth,  the  sum  (p*'c  converges  distributionally 
to  a  function  ip  whose  corresponding  ip*1  is  bounded  and  whose  Fourier  transform  does  not  vanish 
at  0.  Then  the  approximation  order  provided  by  ip  is  (exactly)  k. 


Proof:  By  the  definition  of  ip, 

(3.7)  p(i-)ip=uK, 

where  u  (the  Fourier  transform  of  c)  is  a  27r2Z<i-periodic  tempered  distribution.  Since  ip*'  is  bounded 
by  assumption,  ip  is  continuous  (recall  the  remarks  following  Proposition  2.3).  Since  K  is  continuous 
and  A'(0)  /  0  (and  the  left  side  of  (3.7)  is  continuous),  we  conclude  that  u  coincides,  on  some 
neighborhood  SI  of  the  origin,  with  some  continuous  function.  Therefore 

ip  =  —77-r  on  2ir7Ld  +  Sl\ 0. 

P(f) 

Since  p(D)  is  elliptic,  |p|  >  const|  •  |fc  for  some  const  >  0,  while  const')  •  \k  >  |p|  for  some  const',  by 
the  homogeneity  of  p.  Thus,  for  small  enough  x  and  for  f)  6  2^2Zd\0,  by  the  periodicity  of  u, 


+  0) 
ip{x) 


.  K(x  +  P)p(ix) 
'p(i(x  +  p))K(xy 


<  constjx|fc|/?|“fc. 


Since  k  >  d,  we  see  that  ip  satisfies  the  synthesis  condition  of  order  k,  hence  (by  the  Lower  Bound 
Theorem  2.37)  the  approximation  order  is  at  least  k. 

Now,  choose  fi  G  27rZZd\0  such  that  A'(/3)  ^  0.  Since  p(D)  is  homogeneous  and  elliptic  of  order 
k,  p(ih6)/hk  ss  ikp{9)  ^  0.  Hence,  for  any  9  €  IRd\0, 

,-J(h9  +  P)  k  K(h9  +  p)p(ih0)  ffc  K((d)p{9) 

$(M)  p(i(h9  +  p))I((hO)  ft-o  p(i/3)A(0)  T  U‘ 

Since  ip  does  not  vanish  on  some  neighborhood  of  the  origin,  the  Upper  Bound  Theorem  2.18  applies 
to  yield  that  the  approximation  order  is  at  most  k.  4 
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The  extension  of  the  above  theorem  to  non-homogeneous  elliptic  operators  is  straightforward, 
and  is  omitted  only  because  of  lack  of  examples  of  this  type  in  the  present  literature.  A  discussion 
of  non-elliptic  cases  appears  in  the  box  spline  section. 

As  mentioned  before,  in  case  ip  :=  T <p  decays  only  slowly  at  oo,  the  lower  bounds  on  the 
approximation  order  now  in  the  literature  usually  underestimate  the  actual  approximation  order. 
To  explain  this,  it  is  instructive  to  compare  the  approximation  scheme  suggested  here  in  the  lower 
bound  analysis  with  the  simpler  scheme 

(3.8)  /(&•) «  ip*'f(h-). 

Assume,  for  simplicity,  that  /  e  V,  hence  f(hx)  =  (2n)~d  JiR<i  f(w)ehs(w)  dw.  The  approximation 
scheme  here  is  of  the  form  f(h-)  «  (2 ir)~dip*'  dw.  Taking  into  account  the  fact 

that  the  scheme  here  is  proved  to  be  optimal  (in  terms  of  approximation  orders),  the  optimality 
of  (3.8)  depends  on  the  behaviour  of  the  difference  ip  -  1  around  the  origin,  or  more  precisely,  on 
the  order  of  the  zero  which  ip  -  1  has  at  the  origin.  Indeed,  the  difference  operator  T  is  meant  to 

/S  A 

produce  a  high  order  zero,  but,  since  <p  is  smooth  away  from  the  origin,  a  high  order  zero  of  ip  -  1 
at  the  origin  implies  that  ip  is  globally  smooth,  hence  ip  decays  fast  at  infinity,  which  is  contrary 
to  the  assumption  here. 

Furthermore,  resolving  the  singularity  of  $  at  the  origin  with  the  aid  of  a  trigonometric  poly¬ 
nomial  or  another  smooth  periodic  function  might  be  hard  in  case  this  singularity  is  not  of  integral 
order  (e.g.,  some  log  singularity  or  fractional  singularity).  This  explains  why  in  some  cases  it 
was  impossible  to  remove  the  log  factor  in  the  approximation  order  by  further  differencing  ([Bu3], 
[DJLR]). 


3.3.  Spectral  approximation  orders 

A  J 

The  analysis  of  the  scaling  case  shows  that,  for  this  case,  the  vanishing  of  <p  on  2?rZo  \0  is 
necessary  for  obtaining  positive  approximation  orders.  However,  high  approximation  orders  can 
at  times  be  obtained  from  spaces  spanned  by  basic  functions  whose  Fourier  transform  vanishes 
nowhere  on  IRd,  even  though  the  function  scale  {4>h}h  involves  the  dilates  of  a  single  function. 

Suppose  that,  for  some  function  (p  and  some  neighborhood  II  of  the  origin, 


..Si ±im 
11  $ 


llioofn/*) 


decays  fast  to  0  not  only  as  0  — ►  oo,  but  also  as  A  — ►  0.  In  this  case,  we  may  choose 


cph:=<P(X{h).) 


for  appropriately  selected  decreasing  {A (ft)}/,,  with  the  only  limit  on  the  approximation  order  being 
the  rate  of  decay  of  ,A.  as  A  -*•  0. 
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To  simplify  the  analysis,  we  assume  throughout  the  discussion  that  <f>  satisfies  the  following 
condition 

(3.9)  p(\,p)  :=  \\—  +£  -hLgo(n,x)  <  c(s,e,Sl)e-M-'y/y ,  fi  G  2^2Zd\0,  X  >  0 

for  some  positive  s,  some  sufficiently  small  £,  and  some  O-neighborhood  SI.  We  have  chosen  this 
particular  condition  since  it  is  satisfied  by  functions  whose  Fourier  transform  decays  exponentially 
at  infinity,  as  the  following  proposition  shows.  Other  decay  rates  can  be  treated  along  the  same 
lines. 

Proposition  3.10.  Assume  that  <f>  satisfies  the  condition 

(3.11)  0  <  Ci  <  $(x)|elrl*  <  ci  <  oo,  x  G  IRd, 

for  some  positive  s.  Then  <f>  satisfies  (3.9)  with  the  same  s  and  for  any  e  <  2n,  provided  SI  =  Be/ 2. 
Proof:  It  follows  from  (3.11)  that,  for  any  x  G  Bn,  /?  G  2x7Zd\0,  and  X  >  0, 

!?((;  + WA) |<m,s„-(M'-|.|')A'  <  C01lste-(W-WA., 

«*/ A)  - 

since  min|y|=|x|(|y+^|*-|y|J)  =  (|/3|  —  |x|)J  —  |x|a  >  (|/?|-2|x|)a,  using  the  fact  that  |/?|  >  2ir  >  2|x|. 
Thus,  for  any  e  <  2ir,  (3.9)  holds  with  SI  =  Bc/2  ■  4 

Assuming  (3.9),  we  define  <f>h  <f>(X(h)‘)  for  some  positive  sequence  {A(/i)}/l.  This  implies 
that,  for  (i  G  2jrZZd\0, 


*A(A),,J)=||&fct%  <  ce-(W-.)-/w, 

<t>h  ’ 

and  we  are  led  to  the  following  result: 

Theorem  3.12.  Assume  that  <fi(t-)*'  is  bounded  for  every  t  >  0,  that  <f>  vanishes  nowhere ,  and 
that  (3.9)  holds  for  some  neighborhood  SI  of  the  origin,  some  e  <  n/2,  and  some  s  >  0.  Define 
4>h  :=  <f>(X(hy).  If 

(3.13)  e—/W  <  c/^  with  a  ;=  _  3 £)a, 

then  the  approximation  order  provided  by  {Sh  =  S(<f>h)}k  is  at  least  m,  and  dist < 
c'hm\\f\\'0  -f  o(hk)  for  every  k-admissible  f  in  this  case.  In  particular,  if  X (h)  —  0(hr)  for  some 
positive  r,  then  distoo(/(/i-),  S(^>/,))  =  o(hk )  for  every  k-admissible  f  and  every  k.  Moreover,  if 
X(h)  =  h,  then  distoo(/(/i’)>  $h)  =  0(e~a/k‘)  for  very  smooth  functions  f  (e.g.,  functions  whose 
Fourier  transform  is  a  compactly  supported  measure). 

Proof:  Since  the  argument  here  will  be  used  also  in  the  sequel,  we  prefer  to  provide  the  main 
part  of  the  proof  in  the  form  of  a  separate  lemma: 
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Lemma  3.14.  Let  {4>h}h  be  a  sequence  of  functions ,  with  corresponding  positive  bounded  sequence 
{X(h)}/,,  which  satisfy  the  following  three  conditions: 

(a)  Each  4>h*'  is  bounded; 

(b)  each  4>h  vanishes  nowhere  in  some  h-independent  neighborhood  fl  of  the  origin; 

(c)  with  £  and  s  as  in  Theorem  3.12, 

(3.15)  ifr  J'V-ln)  <  0  6  2>rZZ>. 

<Ph 

Then  {(f>h}h  satisfies  the  conclusions  of  Theorem  3.12. 

Proof  of  the  Lemma:  For  (3  €  2tt2Z<1\0,  let  K(f3)  be  the  open  ball  of  radius  £  <  tt/2  centered 
at  (3  -  2£jj|j.  Then  the  balls  K{(3),  (3  £  2?rEd\0,  while  sharing  the  same  volume  6,  are  pairwise 
disjoint,  and  their  union  is  disjoint  from  i?2*-3«.  Furthermore,  with  F(x)  :=  e~lrl’/A  for  some 
A  >  0,  we  have 

inf  F(x)  = 

zei<(p) 

Therefore,  for  A  <  Ao, 

Y  c-(w-c)*/A<  i/6  y  [  F{t)dt 

PGUTLd\e  P&T'U'l\0J 

„1(rt  <  1/6  /  F(t)dt 

=  c(d)  r  rd-xe-T‘ I x  dr 
«/2ir-3e 

=  cM,s,A0)e-<2,r-3e)*'\ 

Using  (3.16),  we  can  now  derive  conditions  on  {A (h)}/,  which  ensure  that  the  //-synthesis  condition 
holds  with  Uh(x)  =  hm.  Since  <f>h  vanishes  nowhere  on  fi  and  (3.15)  holds,  then,  in  view  of  (3.16),  in 
order  to  check  that  the  //-synthesis  condition  holds,  we  need  only  to  verify  that,  for  a  =  (2x  -  3s)* 
and  some  constant  c  >  0  and  for  all  small  enough  h, 

e-*IKhY  <  chm, 


which  is  (3.13),  and  the  first  claim  of  Theorem  3.12  then  follows  from  Corollary  2.39.  The  second 
claim  follows  .from  the  first,  since  it  is  clear  that  (3.13)  holds  for  any  positive  m  (and  sufficiently 
small  h)  in  case  A (h)  =  hr  for  some  r  >  0. 

Finally,  in  the  case  A(/i)  =  h ,  the  //-synthesis  condition  holds  even  for  Vh{x)  =  e~ath> .  Substi¬ 
tuting  this  into  Corollary  2.39,  and  recalling  that  the  o(hk)  term  there  is  the  first  term  of  (2.31), 
we  obtain  the  desired  result  since,  for  a  band  limited  function,  the  first  term  of  (2.31)  is  0  for  small 
enough  h.  4k 

To  complete  the  proof  of  the  theorem,  we  need  to  show  that  the  assumptions  on  {<f>h}h  made 
in  the  theorem  are  stronger  than  those  assumed  in  the  lemma.  Since  any  <ph  in  the  theorem  is  <£(A>) 
for  some  A,  conditions  (a)  and  (b)  in  the  lemma  follow  respectively  from  the  assumptions  in  the 
theorem  that  <f>{t-)*'  is  bounded  for  any  t  >  0,  and  that  <f>  vanishes  nowhere.  As  to  condition  (c)  in 
the  Lemma,  it  is  implied  by  (3.9),  since  now  4>h  =  c^(-/A(/i)).  4k 
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Proposition  3.10  describes  a  simple  condition  which  implies  (3.9).  A  simple  condition  which 
guarantees  the  boundedness  of  <p(t-)*'  for  every  t  >  0  is  that  |<£(x)|  =  0(|x|-*d+^)  for  some  6  >  0, 
as  x  -»  oo.  Thus  we  conclude  the  following. 

Corollary  3.17.  Assume  that  |<£(x)|  =  0(|x|“^+^)  at  oo  for  some  6  >  0,  and  that  (3.11)  holds 
for  some  s  >  0.  Then  the  requirements  of  Theorem  3.12,  whence  its  conclusions,  hold. 

So  far,  we  have  considered  the  case  when  <f>*'  is  bounded.  However,  the  analysis  here  applies 
also  to  cases  when  <p*'  is  not  bounded,  but  boundedness  can  be  obtained  by  localization,  i.e., 
by  differencing  the  original  cp .  In  this  process  a  certain  precaution  is  required:  if  <p  is  the  non- 
localized  function  and  ip  :=  Tip  is  its  localization  (with  T  involving  only  integer  translates),  then, 
in  general,  it  is  not  desirable  to  define  our  sequence  {<ph}h  by  scaling  ip,  i.e.,  a  definition  of  the  form 
<ph  :=  ip{\(h)-)  should  be  avoided.  The  reason  for  this  is  that  in  this  way  we  scale  also  the  difference 
operator,  hence  obtain  a  difference  operator  that  involves  non-integer  translates  (see  the  box  spline 
discussion,  in  which  the  deteriorating  effect  of  non-integer  translates  on  approximation  orders  is 
detailed;  furthermore,  the  zeros  of  the  Fourier  transform  of  T  then  prevent  us  from  finding  a  domain 
ft  where  none  of  the  <p/,  vanish).  What  can  be  done  is  to  generate  first  a  sequence  <p/,  :=  <p{\{h)-), 
and  then  to  localize  each  <ph  separately  with  the  aid  of  a  difference  operator  T/,  (each  of  which 
employs  only  integer  translates).  Here  is  a  typical  result  in  this  direction. 

Corollary  3.18.  Assume  that,  for  a  given  <p,  <p  coincides  with  some  non-vanishing  function  on 
IRAO,  and  that  l/<p  extends  to  a  continuous  function  on  all  of  JR1*.  Assume  further  that  <p  satisfies 
(3.9)  for  some  s,  e  and  ft.  Let  X(h)  be  decreasing  in  h,  and  let  {T^}^  be  a  set  of  difference 
operators  (each  using  only  integer  translates)  such  that,  with  <pk  Th<P{X(h)-),  the  operators  fa*' 
are  bounded,  and  such  that  the  functions  fa  vanish  nowhere  on  ft.  Then,  all  the  results  stated  in 
Theorem  3.12  hold  with  respect  to  this  {fa}/,. 

Proof:  The  result  follows  from  Lemma  3.14  as  soon  as  we  verify  all  the  conditions  specified  there. 
The  boundedness  of  each  fa*  as  well  as  the  nonvanishing  of  fa  on  ft  are  assumed  here.  It  remains 
therefore  to  consider  the  ratios 

<Ph( •  +  P) 

<Ph 

in  order  to  verify  (3.15).  Let  U/,  be  the  Fourier  transform  of  T/,.  Then  u/,  is  27r2Zd-periodic,  and 
therefore  -  iL^PJMPTl  ^  Thus,  (3.15)  is  implied  by  (3.9),  which  is  assumed  here.  4 

Specific  basic  functions,  for  which  one  is  able  to  obtain  infinite  approximation  orders  are 
discussed  in  [M],  [MN2]  and  [BuDl,2],  and  our  interest  in  this  topic  was  stimulated  by  a  discussion 
with  N.  Dyn.  We  remark  that  Madych  and  Nelson  derived  their  results  in  the  more  general  context 
of  scattered  translates. 


3.4.  Box  splines 

Let  5  6  lR‘lxm  be  of  full  rank  d  and  with  no  0  column.  We  will  also  consider  E  as  the  multiset 
{£  :  £  €  5}  of  its  columns  and  therefore  mean  by  Y  C  5  that  Y  is  a  matrix  obtained  from  E  by 
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omitting  some  columns.  Let  A  :=  {A^}^g=  be  an  arbitrary  set  of  complex  scalars.  The  (exponential) 
box  spline  M  :=  is  defined  via  its  Fourier  transform  as  follows: 


(3.19) 


M(w)  :=  TT 
teBJ  o 


i.e.,  it  is  the  convolution  product  of  the  functionals 

Jo 

In  general,  M  is  a  measure  supported  in  E[0..1]m.  Since  we  assumed  that  rankE  =  d,  the  box 
spline  M  is  a  bounded  function,  of  global  smoothness  fc(E)  -  2,  with 


(3.20)  k  :=  k(E)  :=  m\n{#Y  :  Y  C  E,  rank(E\K)  <  d). 


For  a  generic  choice  of  A  £  Cm,  M  is  a  piecewise-exponential  function  (called  an  exponential 
box  spline,  or  simple  exponential  box  spline);  otherwise,  it  is  a  piecewise-exponential-polynomial 
function.  Piecewise-polynomials  are  obtained  by  the  choice  A  =  0  (polynomial  box  splines,  or  box 
splines).  As  might  be  anticipated  from  their  definition,  box  splines  are  obtained  by  differencing  a 
specific  fundamental  solution  of  the  equation  p(D)f  =  6,  with 

p{x)  :=  H(A«  “(■*), 

and  thus  (cf.  the  discussion  in  the  section  on  the  differencing  case),  polynomial  box  splines  are 
refined  by  scaling,  i.e.,  Af/,  :=  M,  while,  for  general  box  splines,  Mh  is  defined  by 

(3.21)  Mh(w)  :=  TT  t  e{hXc-*  w)t  dt. 

tezJo 

The  point  of  this  refinement  is  that  in  this  way  the  local  structure  is  preserved,  i.e.,  the  pieces  of 
{Mh{-/h)}h  (for  fixed  E  and  A)  all  belong  to  the  same  finite-dimensional  exponential-polynomial 
space.  It  is  important  to  note  that  box  splines  fall  into  the  differencing  case  (as  defined  in  §3.2)  only 
when  E  £  7Ldxm,  since  otherwise  the  difference  operator  used  in  the  localization  of  the  fundamental 
solution  employs  translations  in  non-integer  directions. 

Results  on  approximation  orders  for  box  splines  can  be  found  in  [BD],  [BH],  (DM1],  (E  £ 
2ZdXm,  A  =  0,  see  also  [BHR]),  [Rl],  [DR],  [LJ]  (same  E,  general  A),  [RS]  (general  E,  A  =  0). 
Neither  upper  bounds  nor  lower  bounds  on  the  approximation  order  are  known  for  general  E  and 
A.  In  what  follows,  we  will  derive  upper  bounds  for  the  approximation  order  of  any  box  spline 
(i.e.,  general  E  and  A),  and,  in  case  the  spline  is  smooth  enough,  will  provide  also  matching  lower 
bounds.  Since  the  integral  case  (i.e.,  E  £  TLdxm)  is  the  one  mostly  explored  in  the  literature,  and 
since  our  results  apply  to  this  case  almost  directly,  we  found  it  instructive  to  begin  with  this  special 
case.  Before  doing  that,  we  remark  that  the  operator  sequence  {Mh*'}h  (-  and  A  fixed)  is  always 
uniformly  bounded.  This  follows  from  the  fact  that  supp  A//,  C  E[0 . .  l]m  and  that  the  functions 
{Affc}*<i  are  uniformly  bounded.  The  latter  claim  can  be  verified  as  follows  ([DR]):  let  N  be  the 
box  spline  associated  with  the  same  direction  set  S,  but  with  A  =  0.  From  the  definition  of  the  box 
spline  Mh,  it  follows  that,  as  a  linear  functional, 

Mh:f~  [  ehXtf(Zt)dt ,  f£C(TLld), 

J[ 
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and  thus,  for  h  <  1, 


|M*(/)|  <  cjATfl/l)  <  eA||Ar||  ||/||i.,(s|o..1|~). 

Consequently,  ||M/,||  <  cx||jV||,  all  h  <  1.  For  later  reference,  we  record  this  fact  below. 
Proposition  3.22.  For  a  box  spline  M,  {Mh*'}h<i  is  uniformly  bounded. 


3.5.  Box  splines:  integral  case 

We  assume  here  E  6  7Ldxm  (and  rank(E)  =  d). 

To  start  with,  we  note  that  (3.21)  implies  that  {1/M/,}*  is  uniformly  bounded  in  a  neighbor¬ 
hood  fl  of  the  origin.  This  means  that,  for  the  analysis  of  upper  bounds  on  the  approximation 
order,  we  may  replace  the  quantities  {M*(h0  +  /?)}*  in  the  Upper  Bound  Theorem  2.18  by  the 
ratios 


(3.23) 


M^hd  +  Q)  =Tr  h(\t-it>0) 

Mh(hO)  l|fc(Ac- *•*)-*•/?’ 


/?  e  2tZ zd\o. 


Here  we  have  used  the  fact  that  £  G  TLd  and  /?  £  2x2Zd  implies  that  f  •  f3  £  2v7L,  i.e.,  that 
e-«€-0  _  ^  \ye  now  fix  8  £  ntrf  such  that  A€  -  if  •  0  £  0,  Vf  £  E.  By  the  definition  (3.20)  of  Jfc(E), 
there  exists  Y  C  E  such  that  rank(F)  <  d ,  #(E \Y)  =  k(E)  and  rank(y  U  {£})  =  d,  Vf  £  (E\F). 
Since  rank(y)  <  d  and  Y  C  S  €  7Ldxm,  we  can  find  /?  £  27r2Z<<\0  such  that  YT(i  =  0.  Then,  since 
rank(K  U  (0)  =  d  for  £  6  (E\y),  we  must  have  f  •  f)  ^  0  such  f .  Therefore, 

L-k(S\  Mh(h0  +  (3)  _  r-r  At  -  •  A  _  TT  •  g 

Invoking  the  Upper  Bound  Theorem  2.18  (which  is  applicable  due  to  Proposition  3.22),  we  thus 
obtain  the  following  result,  which  was  first  established  in  [LJ]  by  different  means: 

Theorem  3.24.  let  {£*}/,  be  the  spline  spaces  spanned  by  the  integer  translates  of  the  box 
splines  {M*}*  (resp.)  as  defined  in  (3.21).  Assume  that  E  is  integral.  Then  the  approximation 
order  provided  by  {£/,}/,  does  not  exceed  the  number  k(E)  defined  in  (3.20). 

The  fact  that  the  approximation  order  from  the  box  spline  spaces  {5*}*  is  at  least  k(E)  was 
proved  in  [DR].  The  following  theorem  reproduces  this  result,  but  only  for  sufficiently  smooth  box 
splines. 

Theorem  3.25.  Let  {5*}  be  box  spline  spaces  as  in  Theorem  3.24.  For  every  /?  €  2x7Ld\0,  define 


(3.26) 


J*:={f€S:  le-fl^O}. 


Assume  that 
(3.27) 


e  n 


1 


<  00. 
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Then  the  approximation  order  provided  by  the  spaces  {Sh}h  is  fc(E). 


Proof:  By  Theorem  3.24,  we  only  need  to  show  that  the  approximation  order  is  at  least  k(S). 
Firstly,  note  that  #Lp  >  fc(E)  for  every  /?  €  2n7Ld\0.  Secondly,  all  factors  in  (3.23)  with  £  •  /3  =  0 
equal  1,  hence 

Mh(w  +  /?)  _  |-r  (hAt  -  t‘£  •  w) 


Since  /3  €  2jtZZ<1\0  and  Lp  C  ZZd\0,  it  follows  that  |£  •  /3|  >  2nr  for  every  £  €  Lp.  Thus,  choosing 
t)  :=  min{i[|j  :  £  €  E},  we  conclude  that  for  small  enough  h  (e.g.,  h  <  min$(7r  -  1)/A^)  and  all 


Mh{w  +  fi)  n  hlAfl  +  1£1M 
Mk(w)  \(-m 


<const(A  +  M)#1»  JJ  jjijj 

K  * p 1 


<const(/i  +  |tn|)^s>  JJ  r-i— 
tel* 


Since  we  assume  (3.27),  it  follows  that  {M/J/,  satisfies  the  synthesis  condition  of  order  Jfc(E), 
and  the  desired  result  thus  follows  from  the  Lower  Bound  Theorem  2.37.  4 


Remark.  It  should  be  noted  that  the  condition  (3.27)  in  the  theorem  above  is  active,  namely, 
it  does  exclude  box  splines  of  low  smoothness,  and  in  this  regard  the  theorem  is  weaker  than  the 
original  result  in  [DR].  For  example,  if  d  =  1,  then  Lp  -  E  for  every  /?  €  2jt2Z\0,  and  therefore 
condition  (3.27)  holds  if  and  only  if  m  :=  #S  >  1.  For  a  general  d,  it  is  easy  to  see  that  the 
condition  k(E)  >  d  +  1  implies  (3.27)  but  not  vice  versa.  It  would  be  nice  to  know  whether  the 
arguments  in  the  theorem  can  be  extended  to  box  splines  of  low  smoothness,  particularly  since  a 
similar  gap  appears  below  in  the  non-integer  extension  of  Theorem  3.25. 


3.6.  Box  splines:  non>integral  directions 

In  case  E  €  2Zdxm,  the  condition  £  •  (3  ^  0  (used  in  the  definition  (3.26)  of  Lp)  is  equivalent 
to  £  'P  €  27r2Z\0,  thus  leading  to  the  simple  formula  (3.23)  for  the  ratio  Mh(hQ  +  P)/Mh{hO).  For 
general  E,  we  obtain  such  a  simple  expression  only  for  the  factors  corresponding  to  £  in 

Kp  :=  {£  6  E  :  £•/?<=  2;rE\0}. 

As  it  turns  out,  the  factors  corresponding  to  £  6  Lp\Kp  are  of  no  help  for  the  approximation  order. 
In  particular,  the  approximation  order  is  now  given  by 

(3.28)  *'(E)  :=  min{#A'0  :  €  27r2Zd\0} 

rather  chan  by  the  possibly  larger  &(E)  =  ic\\n{#Lp  :  /?  €  27r5Zd\0}  (cf.,  (3.20)  and  (3.26)).  The 
following  result  has  been  proved  in  [RS]: 
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Result  3.29.  The  approximation  order  provided  by  any  polynomial  box  spline  M  (i.e.,  A  =  0) 
equals  the  integer  k'(E)  defined  in  (3.28). 

Result  3.29  generalizes  the  original  result  from  [BH]  in  which  polynomial  box  splines  with 
integral  direction  set  were  considered.  Both  proofs  are  based  on  the  identification  of  the  polynomials 
in  S(M )  (cf.  Result  3.2).  Furthermore,  the  extensions  of  the  result  of  [BH]  to  arbitrary  A’s  made 
use  of  the  exponential- polynomial  space  in  S  (lower  bounds,  [DR]),  and  ti;e  local  structure  of  the 
box  spline  M  (upper  bounds,  [LJ]).  However,  for  a  non-integral  E  and  nonzero  A,  the  study  of  the 
exponential-polynomials  in  S  or  the  local  structure  of  M  seems  to  fall  short:  the  approximation 
properties  of  the  relevant  exponential-polynomial  space  provide  lower  bounds  on  the  approximation 
order  which,  in  some  cases,  underestimate  the  correct  order,  and  the  local  structure  of  M  provides 
upper  bounds  on  the  approximation  order  which,  usually,  overestimate  the  correct  order. 

In  contrast,  the  analysis  here  of  the  approximation  orders  for  general  box  splines  follows  the 
outline  of  the  analysis,  given  in  the  preceding  section,  for  a  box  spline  with  integer  directions,  with 
some  necessary  modifications.  In  this  way,  the  results  of  the  preceding  section  are  shown  to  hold  for 
a  general  E,  with  fc'(E)  replacing  k( E)  (which  is  a  true  generalization  since  we  have  k'( E)  =  k( S) 
for  E  €  7Ld*m). 

Let  (3  €  2ir7Zd\0.  A  straightforward  computation  then  shows  that,  for  any  8  6  IRd  and 

f  6  s\ 

[  e(hh-H-(P+M))t  dt  — ♦  [  e~^‘0tdt^  0. 

Jo  /l“*°  Jo 

On  the  other  hand,  if  f  6  Kp  and  6  €  satisfies  A  {  -  if  •  0  ^  0,  then  we  get  that,  for  this  f , 

a-1  /'  dt  _  k^iLi 

If  we  now  choose  €  27r2Zrf\0  for  which  #Kp  =  k'(E),  and  choose  6  €  IRrf  such  that  A$  -  if  •  6  £  0 
for  every  i ;  €  Kp,  we  conclude  from  the  above  arguments  that 

I  Mh(h0  +  (3)\to(hk'&). 

Furthermore,  Proposition  3  22  ensures  that  {Mh*')h  is  uniformly  bounded,  and  therefore,  via  the 
Upper  Bound  Theorem  2.18,  we  finally  arrive  at  the  following  generalization  of  Theorem  3.24. 

Theorem  3.30.  For  each  h  >  0,  let  Sh  =  5(Af/,)  be  the  spline  space  spanned  by  the  integer 
translates  of  the  box  spline  Mh  defined  in  (3.21).  Then  the  approximation  order  provided  by  {5/,}/, 
does  not  exceed  k'(E)  (as  given  in  (3.28)). 

In  an  analogous  way,  Theorem  3.25  can  be  generalized  to  non-integral  E  as  follows: 

Theorem  3.31.  Let  {Sh}  be  box  spline  spaces  as  in  Theorem  3.30.  For  every  (3  6  2it7Ld\0,  define 

(3.32)  !/»:={*€  2:  K-/3|^0}. 
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Assume  further  that  3  6  QdXm  and  that 


(3.33) 


e  n 

/3€2nZd\0  ieLf 


1 


<  00. 


Then  the  approximation  order  provided  by  the  spaces  {Shjh  is  k'(3). 

Proof:  Since  we  know,  by  Theorem  3.30,  that  the  approximation  order  is  at  most  k'(3),  we  need 
only  to  prove  that  it  is  at  least  k'(3).  Since  is  uniformly  bounded  in  some  neighborhood  of 

the  origin  (for  small  enough  /i),  it  suffices,  for  an  application  of  the  Lower  Bound  Theorem  2.37,  to 
consider  the  quantities  M^(*  +  /3),  (3  6  27rZd\0  and  h  small  (rather  than  the  ratio  Mh(-  +  (3)/M )»). 
We  now  fix  /3  €  27rZSli\0  and  consider  three  different  cases  of  f  6  3: 

Case  f •  (3  =  0  (i.e.,  f  €  E\L/ j):  In  this  case,  for  sufficiently  small  h  and  for  every  w  6  IRd, 
we  have 

|  /  C(AA< dt\  <  2. 

Jo 

Case  f  •  13  g  2tTL  (i.e,  f  €  Lp\I(p):  Here  we  use  the  estimate 


(3.34) 


|  {  e(h\t-xtiw+0))t  < 

Jo 


3 _ 

|/iA{  -  if  ■  ( w  +  /?)|  ’ 


valid  for  all  w  and  sufficiently  small  h.  Now,  since  E  €  Qdxm,  there  exists  n  €  2Z  such  that 
nE  €  Z£dxm,  and  thus,  since  f  •  ft  ^  0  (and  /?  €  2:rZZrf\0),  we  have  |f  •  /?|  >  /n.  Thus,  (3.34) 

shows  that  for  sufficiently  small  h  and  w, 

Vo  1  If  )®i' 


Case  f  •/?  €  27r2Z\0  (i.e,  f  6  Kp):  In  this  final  case, 


(3.35) 


|  /  e(fc*c dt\  _ 

Jo 


|eu<-i*.w  _  ^ 
|fcA*-if  .(u>  +  /?)|* 


The  denominator  in  the  right  hand  side  of  (3.35)  can  be  estimated  as  in  the  previous  case,  while 
the  numerator,  for  sufficiently  small  h  and  |w|,  can  be  bounded  by  c(h  +  |u>|),  hence  we  finally 
obtain  in  this  case  the  estimate 


if— 

Combining  all  these  estimates,  we  obtain  that,  for  some  /^-independent  neighborhood  Q.  of  the 
origin  and  for  some  /^-independent  constant  c  we  have,  for  small  enough  h  (where  the  “smallness” 
of  h  is  again  /^-independent)  and  all  w  6  il, 


ifi»(»+«i<c(*+M)*''*  n  iAi 

<6 Lf  K 


<  c(h  +  |u7|)*,(s)  PJ 


1 

K-/?r 
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Application  of  the  Lower  Bound  Theorem  2.37,  in  view  of  the  assumption  (3.33),  then  completes 
the  proof.  4 


4.  Examples 


Example  4.1.  Assume  that  <j)  satisfies  the  following  conditions:  (a)  <j>  6  C(lRrf\0);  (b)  for  some 
6  >  0,  |^(x)|  =  0(|a:|~<1~i)  as  x  -+  oo;  (c)  for  some  p  >  0,  |0(x)|  ~  |x|-/1  as  x  -►  0.  Under  various 
additional  conditions  on  <f> ,  (<f>  is  radially  symmetric,  4>  vanishes  nowhere,  and  others;  cf.  [Bu3;pp. 
72-74]),  it  is  shown  in  [Bu3]  that  there  exists  c  :  7Ld  -*■  C  such  that  <f>*'c  converges  absolutely  and 
uniformly  on  compact  sets  to  a  function  x  which  decays  like  at  oo,  and  which  is  a  cardinal 

function ,  i.e.,  satisfies  x(a)  =  ^a,o>  <*  6  2Zd.  [Bu3]  then  proceeds  to  show  that  the  error  in  the 
approximation  scheme 

f(h-)  w  x*'f(h-) 

behaves  like  0(/i-/1)  or  0(h~^\ logh|),  depending  on  the  decay  rate  of  x  at  oo  and  the  type  of  p 
(non-integer,  integer,  even  integer),  and  that  in  some  cases  these  rates  are  sharp  (for  the  above 
approximation  scheme);  cf.  [Bu3;Cor.  5-12].  Ovr  result  in  this  regard  is  the  following: 

Theorem  4.2.  Assume  that  x  is  a  function  for  which  x*'  is  bounded,  and  x(0)  ^  0.  Assume 
further  that  x  can  be  factored, 

X  = 

with  u  a  2n7Zd -periodic  distribution  which  coincides  with  some  bounded  continuous  function  on  a 
neighborhood  fl  of  the  origin,  and  <j>  a  distribution  which  coincides  with  a  continuous  function  on 
(2 irZZd  +  fi)\0.  Assume,  finally,  that  <f>  satisfies  the  following  conditions: 

(4-3)  53  ||0('  +  ^)lk«>(n)  <  oo; 

P  62jr2Z‘<\0 

and 

c\$(x)\  >  |x|~ x  e  fl. 

J  ^ 

Then  the  approximation  order  provided  by  x  is  at  least  p.  Moreover,  if  for  some  6  6  IR  \0,  h^^hd) 
is  bounded  for  small  enough  h,  and  </>(/?)  ^  0  for  some  (3  6  2ir7Ld\0,  then  the  approximation  order 
provided  is  exactly  p. 

Proof:  The  proof  is  very  similar  to  that  of  Theorem  3.6.  By  assumption,  <f>  is  a  well-defined 
function  on  (27rZZd  -f  fi)\0,  and 

»TFf  "1-(n)  =  nfiprf 1 »*-«» s  ci*+«n<..<°>- 

In  view  of  assumption  (4.3),  we  therefore  conclude  that  the  synthesis  condition  of  order  p  holds 
here,  and  consequently,  by  the  Lower  Bound  Theorem  2.37,  the  approximation  order  is  at  least  p. 
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For  the  sake  of  upper  bounds,  we  adopt  the  additional  assumptions  in  the  theorem,  and  note 
that  then,  for  the  same  (3  as  in  the  theorem,  we  have  |  ■■■  1  >  c  >  0,  for  small  enough  h.  On 

4>(h6) 

the  other  hand,  I-  =  | k  I,  and  s‘nce  X  >s  continuous  and  does  not  vanish  at  0, 

\  2  conclude  that  |h-,1x(/i0  +  /?)|  >  ci  >  0.  Application  of  the  Upper  Bound  Theorem  2.18  implies 
that  the  approximation  order  is  at  most  fi.  4 

As  in  the  differencing  case,  the  above  result  shows  that  approximation  orders  provided  by  S(x) 
are  essentially  independent  of  the  localization  process  (i.e.,  the  type  of  periodic  function  u  which 
connects  4>  and  x).  The  only  requirements  on  x  are  that  x(0)  ^  0  and  that  x*'  be  bounded.  For 
this  reason,  we  will  refer  in  what  follows  to  such  approximation  orders  as  “the  approximation  orders 
provided  by  <t>n.  Here  are  some  specific  examples  which  are  covered  by  the  last  theorem.  These 
examples  were  treated  by  several  authors  under  various  restrictions  on  the  underlying  parameter  7 
(see  below)  and  the  parity  d  of  the  spatial  dimension.  A  partial  list  of  references  includes  [J],  [M], 
[MN1,2],  [DJLR],  [Bul-4]  and  [P].  These  references  provide  localization  processes  as  well  as  lower 
bounds  on  the  approximation  order,  the  latter  being  of  the  form  d  -f  A:  or  0(hd+k\\ogh\)  for  some 
k  <  7  whose  value  depends  on  the  quality  of  the  localization  process.  Here  are  the  details. 

(1) :  <j>  =  |  •  |7,  7  £  IR+\22Z+.  Since  $  =  c\  •  |"d“ 7  in  the  complement  of  the  origin,  4>  satisfies 
the  assumptions  of  Theorem  4.2  with  ft  :=  d  +  7,  hence,  indeed,  this  theorem  shows  that,  whatever 
localization  process  is  employed,  the  approximation  order  is  exactly  d  +  7. 

(2) :  <f>  =  |  •  |7  log  |  •  I,  7  €  2IN.  In  this  case,  <f>  coincides  on  lRti\0  with  the  function  c|  •  |“^+^ 
(and  so,  this  case  complements  the  previous  one).  Since  <j>  satisfies  the  assumptions  of  Theorem 
4.2,  with  (i d  +  7,  this  theorem  shows  that  d  +  7  is  the  exact  approximation  order,  regardless  of 
the  localization  process  and  the  parity  of  d. 

(3) :  <fi  :=  (|  •  |2  +  A2)7/2,  7  >  - d ,  7  £  2Z£+.  Here,  <f>  -  c\- 1 -(d+7)/y',  with  c  0,  K  continuous, 
vanishing  nowhere,  and  decaying  exponentially  to  0  at  00.  Thus,  (f>  satisfies  the  assumptions  of 
Theorem  4.2  with  n  :=  d  +  7.  Further  consideration  of  this  </>,  in  which  A  depends  on  h  to  obtain 
spectral  orders,  will  be  given  later. 

(4) :  <f>  —  (|  •  |2  +  A2)7/2log(|  •  |2  +  A2),  7  €  2 7Z+,  d  even.  In  this  case,  4>  admits  a  similar 
expression  to  that  of  the  previous  case,  and  again  we  obtain  approximation  order  d  +  7. 

Example  4.4.  We  discuss  here  several  examples  in  which  spectral  approximation  orders  are  ob¬ 
tained. 

(1)  The  Gaussian  kernel:  <f>  =  /4.  It  was  proved  in  [Bu3]  that  cardinal  interpolation 

using  this  function  does  not  allow  reproduction  of  any  polynomials.  Theorem  3.1  here  shows  that 
further,  since  <f>  =  ce~ vanishes  nowhere,  the  approximation  order  provided  by  cf>  is  0.  However, 
employing  Corollary  3.17,  we  see  that  Theorem  3.12  applies  here  (with  $  =  2).  In  particular,  we 
obtain  the  following  result,  which  corresponds  to  the  choice  A (h)  =  y/h  in  Theorem  3.12. 

Corollary  4.5.  For  the  choice  <j>h  :=  we  have  dist00(/(/i>),  5/,)  <  Ckhk\\f\\'0  +  o(hk )  for 

every  k- admissible  f  and  for  every  k. 

Note  that,  in  terms  of  the  original  basic  functions  {il>h}h  (i.e.,  prior  to  the  scale-up),  the  above 
corollary  suggests  the  choice  ipt ,  =  !k. 
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An  identical  analysis  can  be  made  with  respect  to  other  smooth  functions.  E.g.,  we  can  take 
<j>  —  (|  •  |2  +  l)“(rf+1)/2.  In  this  case,  $  =  ce“l*l,  and  Corollary  3.17,  hence  Theorem  3.12,  apply, 
with  s  =  1,  so  that  we  get,  for  example,  the  following  result. 

Corollary  4.6.  The  conclusion  of  Corollary  4.5  holds  also  for  <j>h  :=  (|  •  |2  +  hr)-(d+1)/2J  for  any 
negative  r. 

Proof:  We  only  need  to  observe,  in  the  application  of  Theorem  3.12,  that,  for  this  (f> ,  <f>( A*)  = 
c(A)(Ma  +  A-*)-W+1)/3.  * 

Example  4.7.  We  continue  with  spectral  orders:  In  the  two  examples  above,  the  operator  <f> *' 
was  bounded.  However,  as  Corollary  3.18  suggests,  our  analysis  also  applies  to  cases  when  <f>*' 
is  not  bounded.  For  instance,  we  can  take  either  ipx  :=  (|  •  |2  +  A-2)7/2,  7  >  -d,  7  £  2ZZ+,  or 
^  =  (I  *  I2  +  A-2r/2log(|  •  |2  +  A’2),  7  €  22Z+,  d  even;  (note  that  these  functions  have  been 
considered  before,  but  in  a  different  context).  In  both  cases  =  c\K(‘/X)\ •  |-d~7  on  !Rd\0,  where 
K  6  Coo(lR‘i\0),  Ji'(x)  ~  e~lIl(l  +  |x|(rf+7“1)/2)  on  all  of  IR/*,  all  derivatives  of  K  of  orders  <  2d+7 
are  in  Li(JR.d ),  and  derivatives  of  I(  of  any  order  (regarded  as  functions  on  IRd\0)  are  rapidly 
decaying  at  00  (all  these  properties  can  be  derived  from  the  known  properties  of  the  modified 
Bessel  functions  [AS],  since  K  =  |  •  |(d+7)/2A(d+-y)/2,  with  I(u  being  the  modified  Bessel  function  of 
third  kind  and  order  u).  Now,  assume  that  <f>h  is  a  localization  of  V»a(/i)j  namely,  <f>h  =  Uh$x(K)  for 
some  27r5Zd-periodic  Uh,  and  <j>h*'  is  bounded.  Since  the  only  singularity  of  V\\(A)  is  at  the  origin, 
we  can  assume  that  uA  does  not  vanish  on  some  punctured  h-independent  neighborhood  fi\0  of 
the  origin  (this,  in  turn,  forces  the  use  of  an  infinite-difference  operator,  except  in  some  special 
cirumstances).  This  ensures  that  $/,  does  not  vanish  on  fl\0,  and  we  further  assume  that  Uh  is 
chosen  such  that  ^(0)  ^  0.  If,  at  this  point,  we  prefer  to  fix  A (h),  i.e.,  if  we  do  not  change  (f>h  with 
ft,  then  we  obtain  a  special  example  of  the  scaling  case,  and  in  such  a  case  the  factor  in  <ph  which 
determines  the  approximation  order  is  |  •  as  can  be  observed  from  Theorem  4.2.  However, 

if  we  change  A (h)  with  h,  as  we  do  in  the  context  of  spectral  orders,  the  dominant  factor  in  <f>h 
becomes  K(-/X(h)).  In  such  a  case,  Corollary  3.18  implies  the  following: 

Corollary  4.8.  Let  {j’\(h)}h  be  as  above,  and  let  <f>h,  h  >  0,  be  a  localization  of  namely, 
<j>h,  =  (h)  for  some  2n7Ld -periodic  function  Uh,  and  <f>h*'  is  bounded.  Assume  further  that 

Uh  vanishes  nowhere  on  some  h-independent  punctured  neighborhood  fi\0  of  the  origin  and  that 
<ph( 0)  0.  Then  all  results  stated  in  Theorem  3.12  hold  with  respect  to  this  {4>h}h  and  with  s  —  1. 

Proof:  We  wish  to  apply  Corollary  3.18,  hence  need  to  verify  that  all  the  conditions  required 
there  hold  in  our  case.  By  assumption,  the  operators  {<f>h*'}h  are  bounded.  Also,  the  various 

A  A  1 

assumptions  on  the  zeros  of  Uh  and  <f>h,  together  with  the  fact  that  <f>  vanishes  nowhere  (on  1R  \0), 
imply  that  (f>h  vanishes  nowhere  on  ft.  Thus,  in  order  to  apply  Corollary  3.18,  it  remains  to  verify 
that  <f>  :=  (|  •  |2  +  I)7/2  (or  4>  :=  (|  •  |2  +  I)7/2  log(|  •  |2  +  1))  satisfies  (3.9).  For  this,  we  note  that  the 
fact  that  K  ~  e”W(l  +  |  •  |(d+-y— l)/2 )  impijes  that 


&(»  +  /?)/A) 

£(x/A) 


<  const 


e-\(x+p)/x\ 

e-l*/*! 
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for  sufficiently  small  x  and  A  and  for  all  (3  £  27r5Zd\0.  This  implies,  as  in  Proposition  3.10,  that  $ 
satisfies  (3.9)  with  s  =  1.  4 

Some  improvements  of  the  above  result  are  available.  First,  we  have  neglected  the  “positive” 
role  of  the  factor  |  •  In  fact,  let  <ph  :=  ipxth)  as  before.  Then,  for  x  6  Bc  and  /?  €  27r2Zd\0 

(and,  say,  with  X(h)  <  1),  the  estimate 


1  +  (1Wi)(d+7~1)/2 

|x  +  (1  +  |x/A(/i)|)(d+-y-i)/2) 


<  const  A (h)^^"1)/2  |/j|-(^+i)/2 


holds  (for  the  case  d  +  7  >  1;  otherwise,  the  factor  A(h)"(ti+7''1^2  can  be  removed  from  the  above 
bound  and  the  subsequent  analysis  becomes  simpler).  Thus,  for  the  case  d  +  7  >  1,  we  get  that 

<j>h(x  ~h  P)  _ 

Mx)  ~  K(x/X(h))\x  +  /3\^ 

<  const  ^  e-^W>Wl(1  +  (^l).^-.,/.) 

_  const  |i|  C_|l/A(fc)|  |x  +  ^d+7  (]_  +  (|x/A(h)|)(d+7-i)/2) 

<  const  |x!d+7  g-d^l -*')/W/X(k)<d+*>-1V2. 


We  can  estimate  the  sum  J^p&2ir7id\o  e'd/3!-2*)/^)  as  jn  the  proof  of  Lemma  3.14  to  obtain 
that,  for  a  :=  2tt  —  6e, 


\\hLUQ. 


LM)  <  const  e-a/xw/X (/i)(d+7~1)/2. 


By  changing  a  if  needed,  and  assuming  that  A (h)  is  small  enough,  we  can  replace  e— “/^(*)/ A(/i)^d+n'-1^/2 
by  which  means  the  i/-synthesis  condition  holds  here  with 


Substituting  these  {vh}h  into  (2.34)  and  combining  (2.34)  with  (2.31),  we  get  the  following: 

Theorem  4.9.  Let  {<f>h}h  be  as  in  Corollary  4.8.  If  the  sequence  {A (h)}/,  satisfies  e~alx^h'1  <  chm 
for  some  non-negative  m  and  some  a  <  2 x,  then  for  every  k-admissible  function  f,  and  every 
j  <  min{d+  7,fc}, 

(4.10)  dist oo(f(h-),S(<t>h))  <  const  hm+i/||'  +  o(hk). 

In  particular,  if  k  >  d  +  7,  then 

dist oo(f(h-),S(<t>h))  <  const  hm+d+,y\\f\\'d+1  +  o(hk). 

Here  the  constant  is  independent  of  f  and  h. 

The  last  theorem  implicitly  suggests  the  optimal  choice  of  { A(/i)}/l:  if  /  is  fc-admissible,  then 
one  should  choose  m  to  be  any  number  larger  than  k  -  min{d  +  7,fc},  for  that  would  make  the  first 
term  in  (4.10)  go  to  0  faster  than  the  second  one. 
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A  second  remark  here  concerns  the  type  of  difference  operators  that  can  be  used  in  the  lo¬ 
calization  process  leading  to  the  sequence  { <ph]h •  As  previously  explained,  we  cannot  scale  Ti  to 
obtain  T/,,  namely  (assuming  without  loss  that  A(l)  =  1),  we  cannot  take  Uh  :=  uj(-/A(/i)).  On 
the  other  hand,  in  the  present  context,  the  singularity  of  <p  at  the  origin  is  homogeneous,  i.e.,  the 

A 

product  of  <p  by  a  homogeneous  function  G  (viz.  |  •  |d+7)  is  a  continuous  function  which  does  not 
vanish  at  0  (nor  anywhere  else,  for  that  matter).  This  suggests  that  the  same  difference  operator 
that  is  used  to  localize  <p ,  can  also  be  used  to  localize  any  scale  of  <p.  Here  is  a  sample  statement  in 
this  direction: 

Proposition  4.11.  Let  <f>  be  as  in  Corollary  4.8,  and  assume  that  d  +  7  >  1.  Let  u  be  a  smooth 
2ir2Zd -periodic  function  such  that  |  •  | ■'('*+'*) u  has  bounded  derivatives  in  a  neighborhood  of  the 
origin  up  to  order  d+1  inclusive.  Assume  that  ip  is  a  function  which  satisfies  ip  =  u<p{-/ A)  for  some 
positive  A.  Then  if)*'  is  bounded. 

Proof:  Recall  that  <p(-/ A)  =  c/v  (*/A)|  •  where  all  derivatives  of  K,  hence  of  K(-/ A),  up  to 

any  order  <  2d +7,  are  in  L\.  Combining  this  with  the  present  assumption  on  u  (and  using  the  fact 
that  2d+7  >  d+ 1  here),  we  conclude  that  the  derivatives  of  if)  =  cuK(-/X)\  •  |“(d+7)  of  orders  up  to 
d+1  are  integrable  around  the  origin.  Further,  away  from  the  origin,  all  derivatives  of  K,  hence  of 
A'(*/A),  are  rapidly  decaying,  while  the  derivatives  of  u  are  bounded  (by  virtue  of  the  smoothness 
and  periodicity  of  this  function)  and  the  derivatives  of  |  •  |-(d+7)  are  of  (at  most)  polynomial  growth. 
This  proves  that  all  derivatives  of  ip  up  to  order  d+1  are  in  /^(IR**).  Consequently,  ip  =  o(|-|“(d+1l) 
at  00,  which  implies  that  ip*'  is  bounded,  as  claimed.  4 
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